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CHAPTER  1 
INTRODUCTION 

Due  to  the  high  strength-to-weight  and  s tif f ness-to- 
weight  ratios,  composite  materials  are  ideal  for  weight- 
sensitive  structures  such  as  aircraft,  spacecraft,  and 
automotive  vehicles.  In  recent  years,  with  the  advent  of 
jet  propulsion,  particularly  with  the  current  increased 
interest  in  short  take-off  and  landing  aircraft,  it  has 
become  necessary  to  pay  increasing  attention  to  the 
higher  frequency  motions  of  such  structures.  These 
motions  depend  strongly  on  the  structure’s  damping  or 
capability  for  dissipation  of  vibratory  energy  [1].  In 
addition,  a  new  type  of  excitation  has  become  more  preva¬ 
lent,  random  excitation  either  of  mechanical  or  acoustical 
origin  [2].  For  example,  jet  engine  exhaust  generally 
contains  a  noise  spectrum  wide  enough  to  excite  most  of 
the  natural  frequencies  encountered  in  aircraft  structures 
[3].  The  natural  resonance  phenomena,  so  produced,  can  be 
very  destructive.  Since  near  resonant  conditions  can  no 
longer  be  avoided  in  many  types  of  structures,  the  maximi¬ 
zation  of  damping  within  a  structural  system  provides  a 
most  useful  concept  in  controlling  resonance  [4]. 
Unfortunately,  as  will  be  shown  in  later  chapters,  high 
damping  is  mostly  coupled  with  low  stiffness,  and  high 
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stiffne5is  is  mostly  coupled  with  low  damping.  Therefore, 
the  optimization  of  damping  and  of  the  stiff ness-to-weight 
ratio  is  a  practical  idea  in  designing  a  proper  composite 
material  to  be  used  in  aircraft  and  space  vehicles. 


1.1  Literature  Survey 


Most  results  from  a  series  of  researches  on  damping 
beginning  in  the  1920's  [5]  indicated  that  damping  is  a 
material  property. 

Kimball  and  Lovell  [5]  experimentally  showed  that 
for  stress  cycles  of  frequency  of  from  two  to  three  a 
minute,  up  to  fifty  a  second  the  frictional  loss  (a  kind  of 
energy  loss)  is  independent  of  the  frequency  but  is  depen¬ 
dent  on  the  amplitude  of  strain  of  the  cycles  for  eighteen 
different  solids,  including  several  metals,  glass,  cellu¬ 
loid,  rubber  and  maple  wood,  when  strain  was  below  the 
elastic  limit.  Crandall  [6]  pointed  out  that  the  values  for 
material  damping  which  will  be  introduced  in  Section  2.3.3 
encountered  in  practice  ranged  from  about  0.00001  to  0.2; 
however,  Lazan  [4]  pointed  out  the  material  damping  ranged 
from  0.001  to  0.1.  And  the  material  damping  depends  on 
both  the  amplitude  and  frequency  of  the  oscillation.  If, 
however,  the  system  is  completely  linear,  then  damping  is 
independent  of  amplitude  [6]. 

Recently,  the  composite  materials  got  more  attention 
in  industrial  application.  Lazan  t4]  gave  a  detailed 
review  on  material  damping  of  materials  and  material 
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composites.  Kume,  Hashimoto  and  Maeda  [7]  used  the  damping- 
stress  function,  derived  by  Lazan  [8],  to  calculate  the 
material  damping  of  cantilever  beams.  They  found  that  low 
order  modes  of  a  cantilever  beam  with  equal  maximum  stress 
amplitude  gave  almost  the  satiie  material  damping,  theoreti¬ 
cally,  and  experimental  results  have  the  same  order  of 
magnitude  as  the  theoretical  results  when  the  maximum  stress 
amplitude  is  less  than  a  certain  value.  Schultz  and  Tsai 
(9]  indicated  that  unidirectional  glass  fiber  reinforced 
composites  exhibit  anisotropic,  linear  viscoelastic 
behavior  when  those  undergoing  small  oscillation  and  that 
damping  increases  in  magnitude  with  change  in  fiber  orien¬ 
tation  angle  with  respect  to  loading  direction  in  the 
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used  the  laminated  plate  theory  and  two-dimensional  energy 
approach  to  predict  the  flexural  damping  of  laminated 
composites.  In  thtir  work,  damping  coefficient  was  deter¬ 
mined  by  free-free  flexural  modes  of  vibration  (12].  Siu 
and  Bert  [13]  discussed  the  vibration  of  composite  plates 
having  material  damping.  Suarez,  Gibson,  and  Deobald  [14] 
observed  the  dependence  of  damping  on  frequency  of  fiber 
reinforced  epoxy  or  polyester.  Gibson  and  Plunkett  [15] 
found  that  for  small  strain,  damping  and  stiffness  are 
independent  of  amplitude  of  strains,  but,  once  the  thres¬ 
hold  strain  is  exceeded  (i.e.,  failure  starts),  the 
resulting  increase  in  damping  is  much  more  significant 
than  the  corresponding  reduction  in  stiffness.  Similar 
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results  were  also  observed  by  Tauchert  and  Hsu  [16).  Bert 
and  Clary  [17,  181  gave  a  complete  review  on  measurement 
and  analysis  of  damping  and  dynamic  stiffness  for  compo¬ 
sites.  The  first  paper  to  optimize  the  damping  of  the 
structure  was  perhaps  that  of  Plunkett  and  Lee  [19].  The 
damping  of  a  beam  is  improved  by  introducing  thin  con¬ 
strained  visco-elastic  layers  on  the  top  and  bottom  of  the 
structure.  These  viscoelastic  layers  are  then  stiffened 
by  properly  designed  constraining  layers. 

Cox  [20]  discussed  the  stress  distribution  in  fibrous 
materials.  Cox's  shear  lag  stress  analysis  was  later  on 
used  to  analyze  the  stress  distribution  of  short-fiber 
coiviposites,  as  in  studies  [21,  22].  Photoelasticity  [23, 
24]  and  finite  element  methods  [25,  26]  were  used  to 
investigate  the  stress  concentration  in  the  matrix  around 
fiber  tips  of  short-fiber  composites.  Strength  of  short- 
fiber  composites  was  analyzed  in  several  studies  [27,  28]. 
Analysis  of  complex  moduli  for  such  kind  of  material  was 
presented  in  studies  [29,  30,  31] .  High  damping  of  short- 
fiber  composites  was  analytically  and  experimentally 
observed  in  references  [31,  32].  Material  damping  of 
randomly  oriented  and  unidirectional  laminar  short-fiber 
composites  has  been  discussed  by  Sun,  Wu,  Chaturvedi,  and 
Gibson  [33,  34]. 


1.2  Scope  of  This  Study 

The  objectives  of  this  study  are  to  analyze  the  mate¬ 
rial  damping  and  to  optimize  the  specific  stiffness  (the 
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ratio  of  the  stiffness  to  the  density)  and  material 
damping  of  continous  and/or  discontinuous  fibers  rein¬ 
forced  laminated  composite  structure  elements.  The  work 
involved  in  this  research  is  briefly  introduced  as 
follows: 

A.  To  develop  a  short-fiber  composite  model  to  determine 
the  moduli  of  short-fiber  composites. 

B.  To  analyze  the  stiffness  and  material  damping  of 
unidirectional  laminar  fiber  composites/  randomly- 
oriented  fiber  composites,  and  certain  kinds  of 
laminated  fiber  composites  through  classical  lami¬ 
nated  plate  theory  approach. 

C.  To  analyze  the  material  damping  of  laminated  fiber 
composites  by  an  energy  approach,  where  a  three- 
dimensional  displacement  finite  element  method  is 
used. 

D.  To  optimize  material  damping  and  the  specific  stiff¬ 
ness  of  laminated  composite  plates. 

1.3  Material  Constants  and  Ranges  of  Design  Parameters 
In  this  study,  four  different  kinds  of  widely-used 
fiber  composites  (i.e.,  glass  epoxy,  Kevlar  epoxy, 
graphite  epoxy,  and  boron  epoxy)  are  involved;  and  much 
interest  is  concentrated  on  graphite  epoxy  and  Kevlar 
epoxy,  because,  generally  speaking,  the  former  has  higher 
stiffness,  while  the  latter  has  higher  damping. 
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In  order  to  compare  with  the  experimental  results  of 
Suarez,  etc.  [35,  36],  the  material  constants  used  in  this 
study  are  the  same  as  those  of  the  experimental  specimens. 
Some  material  constants,  which  are  not  given  in  those 
experimental  data,  are  obtained  from  reference  [37]. 

Unless  specially  specified,  the  material  constants  used  in 
this  study  are  given  in  Table  1.1. 

The  length  of  fiber  is  one  of  the  characteristic 
parameters  of  short-fiber  composites.  Due  to  the 
existence  of  a  critical  fiber  length,  s^^  (i.e.  the  minimum 
fiber  length  in  which  the  ultimate  strength  can  be 
achieved  [38]),  there  is  a  minimum  value  for  fiber  length, 


where  Xy  is  the  matrix  yield  stress  in  shear,  and  d  is  the 
fiber  diameter.  The  lowest  of  those  four  fiber  com¬ 
posites  is  2750  MPa,  as  given  in  reference  [37],  for 
graphite  T-300.  And  the  matrix  yield  shear  stress  is, 
according  to  the  manufacturer's  test  results,  97  MPa  of 
AS4/3501  graphite-epoxy  tape.  Therefore,  the  fiber  aspect 
ratio,  the  ratio  of  fiber  length  to  its  diameter,  should 
be 


s/d  2  14.2 


1.2 
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In  this  study,  the  fiber  aspect  ratio  is  chosen  to  be 
between  25  and  lOOOO;  while  the  fiber  volume  fraction 
is  chosen  to  be  0.65  for  the  most  cases  or  0.5  for 
randomly  oriented  fiber  composites. 

Table  J.  .1:  Material  Properties  Data  of  the  Matrix  and 
the  Fibers 


Constants 

Matrix 

Fibers 

Epoxy 

Glass 

Kevlar 

Graphite 

Boron 

(Gpa) 

3.94 

72-4 

99.8 

175.8 

381.9 

E2  (Gpa) 

HHI 

13.8 

6.9 

13.8 

35.0 

G2^2  (Gpa) 

1.465 

27.6 

13.8 

27.6 

70.0 

CM 

> 

0.345 

0.22 

0.376 

0.16 

0.21 

k  (Gpa) 

4.236 

— 

— 

— 

— 

^1 

0.015 

0.0015 

O.Oll 

0.0015 

0.0019 

p(kq/m^) 

1220. 

2539. 

1479. 

1760. 

2481. 

^1 

— 

2 

2 

2 

2 

^2 

— 

1 

1 

1 

1 

] 

I 


where  and  ^2  constants  used  in  Helpin-Tsai 
Equation,  ni  is  the  longitudinal  damping,  k  is  the  bulk 
modulus,  and  p  is  the  density. 


CHAPTER  2 
DAMPING 

2.1  Definition  of  Damping 

The  process  by  which  vibration  steadily  diminishes  in 
amplitude  is  called  damping.  In  many  ways,  the  assumption 
that  systems  possess  no  damping  is  a  mathematical  conven¬ 
ience,  rather  than  a  reflection  of  physical  evidence.  In 
fact,  if  a  system  is  set  in  motion  and  allowed  to  vibrate 
freely,  the  vibration  will  eventually  die  out;  the  rate  of 
decay  depends  on  the  amount  of  damping.  This  reduction  in 
vibrating  amplitude  occurs  because  the  energy  of  the 
vibrating  system  is  dissipated  as  friction  or  heat  or  is 
transmitted  as  sound  [39].  And  this  is  why  damping  is  also 
interpreted  as  any  phenomenon  within  the  body  of  the 
material  where  energy  is  dissipated  [40], 

The  concept  of  an  undamped  system  serves  not  only  a 
useful  purpose  in  analysis,  but  can  also  be  justified  in 
certain  circumstances.  For  example,  if  the  damping  is 
small  and  one  is  interested  in  the  free  vibration  of  a 
system  over  a  short  interval  time,  there  may  not  be 
sufficient  time  for  the  effect  of  damping  to  become 
noticeable.  Similarly,  for  small  damping  one  may  not  be 
able  to  notice  the  effect  of  damping  in  the  case  of  a 
system  with  harmonic  excitation,  provided  the  driving 
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frequency  is  not  in  the  neighborhood  of  any  of  the  natural 
frequencies  of  the  system  (41J.  On  the  other  hand,  damp¬ 
ing  of  a  given  system  should  be  considered  if  this  system 
is  subjected  to  vibration  near  its  resonant  frequencies 
because  damping  has  a  large  influence  on  the  amplitude  in 
the  frequency  rc«jion  near  resonance  [42]. 

2.2  Damping  Mechanism 

There  are  many  mathematical  models  representing  damp¬ 
ing.  The  mechanism  of  damping  can  take  any  of  these 
forms  and  often  more  than  one  form  may  be  present  at  a 
time.  Therefore,  in  order  to  analyze  or  predict  the 
damping  of  a  given  system,  one  ideally  should  take  into 
account  ail  possible  damping  mechanisms;  fortunately,  in 
most  practical  cases,  one  or  two  mechanisms  predominate  so 
that  one  may  neglect  the  effect  of  all  others.  Three 
widely  used  mathematical  models  for  damping  are  introduced 
below; 


2.2 
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.1  Viscous  Damping 

The  viscous  damping  force  is  defined  as 


I 

I 


I 

I 


I 

I 


-CX 


2.1 


where  the  constant  C  (of  dimension  force  per  unit  velocity) 
is  called  the  coefficient  of  viscous  damping.  This  type  of 
damping  occurs  in  lubricated  sliding  surfaces,  dash-pots, 
hydrolic  shock-absorbers  [43].  The  minus  sign  indicates 


that  this  damping  force  is  always  opposite  to  the  direc¬ 
tion  of  the  motion.  The  work  done  by  damping  force, 
namely,  dissipated  energy  during  one  cycle  of 
harmonic  motion,  x  =  A^sincot,  wil  be 


2ti 

f  dx^  2 

<“D>cyc  “]f<J  C  (gj)  dt  =  bCAqU,  2.2 

where  w  is  the  circular  frequency  in  radians  per  unit  time. 
Apparently,  dissipated  energy  due  to  viscous  damping  in  a 
cyclic  motion  is  proportional  to  the  frequency  and  to  thr 
square  of  the  amplitude  of  motion. 

2.2.2  Dry  or  Coulomb  Dcimpinq 

This  type  of  damping  occurs  in  the  sliding  of  dry 
surfaces.  The  damping  force  during  motion  is  constant  and 
is  given,  according  to  Coulomb's  law,  by 

Fjj  =  -fN  2.3 

where  N  is  the  normal  component  of  the  force  upon  the  sur¬ 
face  of  contact,  and  f  is  the  coefficient  of  dry  friction. 
Damping  induced  by  the  joints  is  mainly  because  of  dry 
damping.  And  it  is  known  that  damping  of  built-up 
structures  (i.e.  structure  made  by  joining  together  skins, 
strings,  frames,  etc.)  could  further  be  caused  by  the 
effect  of  the  joint  [1]. 


■tv 
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2.2.3  Material  Dampim 

This  kind  of  damping  is  also  referred  to  as  internal, 
hysteresis  or  structual  damping.  It  is  caused  by  the 
internal  friction,  the  viscoelastic  behavior  of  the 
material,  and  the  interfacial  slip  in  the  material  itself. 

It  is  well  known  that  an  elastic  body  which  is 
repeatedly  stressed  becomes  hot.  If  an  elastic  body  is 
subjected  to  forced  oscillation,  a  positive  work  of  the 
exciting  force  must  be  spent  to  keep  the  amplitude  of  the 
oscillations  constant  in  time.  The  reasi_  i  for  the  heating 
of  the  body  and  the  expenditure  of  external  work  is  the 
internal  friction  of  the  material.  Although  this 
explanation  can  be  easily  accepted  in  a  qualitative  way, 


it  is  more  difficult  to  translate  the  problem  into 
mathematical  terms. 

Three  principal  hypotheses  have  been  proposed  to 
explain  the  phenomenon  of  internal  friction,  i.e.,  the  vis¬ 
cous  theory,  the  hereditary  theory,  and  the  hysteresis  loop. 

Viscous  cheory.  The  viscous  theory  assumes  that  in 
solid  bodies,  there  exist  some  viscous  actions  which  can 
be  compared  to  the  viscosity  of  fluids.  These  viscosity 
effects  are  assumed  to  be  proportional  to  the  first 
time  derivative  of  strain.  The  coefficient  of  the 
proportionality  (constant  for  each  material  at  constant 
temperature)  is  called  the  coefficient  of  viscosity. 

Based  on  this  assumption,  many  mathematical  models  (scch 
as  Maxwell  model,  Kelvin-Voigt  model,  and  three-parameter 
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models  (44,  45])  have  been  introduced  to  represent  dif¬ 
ferent  materials.  For  the  case  of  Kelvin-Voigt  model 
under  normal  deformation,  the  relationship  between  stress 
a  and  strain  c  is  expressed  by 

o  =  ^  2.4 

where  E  is  Young's  modulus  and  C  is  the  coefficient  of 
viscosity  of  the  material. 

Hereditary  theory.  The  hereditary  theory  attributes 
the  dissipation  of  energy  due  to  material  damping  to  the 
elastic  delay  by  which  the  deformation  lags  behind  the 
applied  force  (43],  According  to  this  theory,  the  defor¬ 
mation  at  a  given  instant,  Instead  of  depending  only  on  the 
actual  applied  stress  at  that  time  as  it  would  if  the 
materials  followed  Hooke's  law,  depends  on  all  the 
stresses  previously  applied  to  the  elastic  body.  The 
stress-strain  relationship  is  given  by 

rt 

o  =  Ee  +  J  (Kt,  t)  e(T)  di  2.5 

where  t  is  the  actural  time  and  x  is  an  instant  of  time 
between  t  =  and  t  =  t.  The  function  4>(t,  x)  is 
called  the  hereditary  kernel  or  memory  function. 

Hysteresis  loop.  For  a  material  under  a  cyclic 
loading,  the  stress-strain  curve  is  a  closed  curve  which 
is  called  the  hysteresis  loop.  The  physical  meaning  of 


model,  the  damping  force  is  assumed  to  be  proportional  to 
velocity  and  inversely  proportional  to  frequency,  i.e. 


X 

cl  0) 


If  an  external  force  is  applied  just  enough  to 
balance  the  damping  force  and  to  maintain  a  simple 
harmonic  motion,  x  =  A^sinut,  then 


=  Kx  t  S  i 


Where  Kx  represents  the  elastic  force  of  the  system;  for 


Equation  2.7,  and  the  plot  of  external  force  F^  as  a 
function  of  displacement  is  a  skewed  ellipse,  as  in 
Figure  2.1. 


)  =  1 

Aq  hA^ 
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Figure  2.1:  Sketch  of  Hysteresis  Loop 


The  work  done  by  damping  force  (dissipated  energy)  in  one 
cycle  (Up) eye  is 


2Tr 


(Ud) 


eye 


dx  = 

J  O 


h 

5<*> 


dt  =  TT  h  A, 


2.9 


Thus,  the  energy  dissipated  in  one  cycle  is  proportional 
only  to  the  square  of  the  amplitude.  This  expression 
agrees  with  the  results  of  experiments  of  Kimball  and 
Lovell  [5]  which  indicate  that  for  a  large  variety  of 
materials  such  as  metals,  glass,  rubber  and  maple  wood, 
subject  to  cyclic  stress  such  that  the  strains  remain 
below  the  elastic  limit,  the  internal  friction  is  entirely 


dependent  on  the  rate  of  strain.  Equation  2.9  also  agrees 
with  Lazan's  notes  [46]  on  dissipated  energy.  For  example, 
at  low  amplitudes  of  stress,  the  dissipated  energy  is 
proportional  to  the  square  of  the  stress  amplitude,  and  the 
hysterestic  loop  is  elliptical  in  form. 

Unlike  homogeneous  materials,  fiber  reinforced  mate¬ 
rials  have  interfaces  between  matrix  and  fiber.  When  a 
fiber  reinforced  composite  is  subject  to  a  tensile  strain 
cycle,  the  high  shear  stress  may  cause  the  fiber  matrix 
interface  to  fail  so  that  energy  is  dissipated  by 
friction  as  the  matrix  slides  over  the  fiuers  {47].  Damp¬ 
ing  is  then  increased  due  to  Interfacial  slips  between 
matrix  and  fiber.  In  this  study,  perfect  bonding  between 
matrix  and  fiber  is  assumed?  conseqeuntly,  interfacial  slip 
is  not  considered. 

2.3  Types  of  Damping  Representations 

Many  different  disciplines  have  been  concerned  with 
damping  measurements,  and  this  has  further  complicated 
nomenclature.  Confusion  has  been  caused  not  only  by  the 
large  variety  of  damping  units  used,  but  also  by  the  lack 
of  unique  definition  for  many  well-accepted  units.  It  is, 
therefore,  desirable  to  review  the  various  damping  units 
currently  used  and  to  indicate  relationships  between  them. 
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2,3.1  Damping  Ratio  (c) 

Figure  2.2  shows  a  single  degree  of  freedom  system 
with  viscous  damping,  excited  by  force  F(t). 


Figure  2.2;  Sketch  of  Viscous  Damping  Model 
Its  differential  equation  of  motion  is  found  to  be 

MX  +  C*x  +  Kx  =  F(t)  2.10 


If  F(t)=0,  one  has  the  homogeneous  differential  equation 
whose  solution  corresponds  physically  to  that  of  free- 
damped  vibration.  The  general  solution  to  this  homogeneous 
equation  is 


X 


-  (c/2m) t 
e 


/ 

(c^e  + 


-bt 
C2e  ) 


2.11 


where 


b 


[  (c/2m) k/mj 


2.12 
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Cl  and  C2  are  constants  to  be  determined  by  initial 
conditions.  In  order  to  have  oscillation,  one  will  expect 
to  have 


2 

(c/2m)  <  k/m  2.13 

Apparently,  there  exists  a  critical  value  c^  for  c, 

2 

when  (c/2m)  equals  k/m.  Damping  ratio  (38],  Q,  is 
defined  as 

r  =  c  2.14 

<=c 


where 


•j 

c^  =  (4mk)  2.15 

2.3.2  Logarithmic  Decrement  (6) 

A  convenient  way  to  determine  the  amount  of  damping 
present  in  a  system  is  to  measure  the  rate  of  decay  of 
free  oscillations.  Logarithmic  decrement  [38]  is  defined 
as  the  natural  logarithm  of  the  ratio  of  any  two 
successive  amplitudes,  as  in  Figure  2.3,  of  a  free 


vibration. 
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Figure  2.3:  Sketch  of  Logarithmic  Decrement 

2.3.3  Loss  Tangent  (tan  t) 

It  is  well  known  that  polymer  behaves  as  viscoelastic 
material,  i.e.,  combining  two  material  properties,  one  of 
which  is  perfectly  elastic,  while  the  second  is  viscous 
fluid  [43].  Let  such  a  viscoelastic  material  be  subject 
to  a  sinusoidal  stress  experiment  at  frequency  u  such 
that  the  period  2ti/u  of  oscillation  is  sufficiently  large 
as  compared  tc  the  transit  time  of  elastic  waves  through 
the  specimen  that  stress  and  strain  can  be  considered 
uniform  throughout  the  test  section.  Under  these  condi¬ 
tions,  the  response  to  a  steady-state  sinusoidal  stress  a 
is  a  steady-state  sinusoidal  strain  e  at  the  same 
frequency  [44],  out  of  phase  by  the  angle  <S>,  e.g. 

o  »  Oq  sintot 


2.17 
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e  =  £q  sin((jjt  -  2.18 

Both  the  response  amplitude  and  the  phase-shift  (or  phase 
angle)  <t>  are  frequency-dependent,  but  in  the  linear  range 
Eq  is  proportional  to  a^.  The  phase  relationships  are 
conveniently  shown  in  the  rotating-vector  representation  of 
simple  harmonic  motion,  as  in  Figi  e  2.4. 


Figure  2.4:  Rotating-vector  Representation  of 
Harmonic  Motion 

The  rotating  vector  OB  of  magnitude  E'co  lags  behind 
the  stress  OA  by  4'  radians.  Stress  OA  may  be  resolved 
into  two  components,  E'Gq  in  phase  with  strain  and 
■n/2  radians  out  of  phase  with  strain,  as  in  Figure  2.4. 

Here  E'  is  the  storage  modulus  and  E"  is  the  loss  modulus. 
The  loss  tangent  tan  4  is  defined  as 

tan  4  =  E"/E'  2.19 

Some  authors  call  tan  4  the  loss  coefficient.  The  ratio 
E"/E'  is  a  measure  of  the  ratio  energy  loss  to  energy 
stored.  For  viscoelastic  material,  the  moduli  are  often 
expressed  in  terms  of  a  complex  number,  called  complex 
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modulus.  In  this  study,  the  superscript  *  is  used  to 
indicate  the  complex  modulus,  for  example; 

E*  =  E'  +  iE"  2.20 


where  i  is  /T. 

2.3.4  Specific  Dcimpinq  Capacity  for  Cyclic  Loading  (4',-) 
The  physical  meaning  of  the  hysteresis  loop  of 
Section  2.2 .3  is  considered  here.  Since  materials  do  not 
behave  in  a  perfectly  elastic  manner  even  at  very  low 
stress  [46],  inelasticity  is  always  present  under  all 
types  of  loading,  although  in  many  cases  extremely  precise 
measurements  are  necessary  to  detect  it.  Under  a  cyclic 
loading  condition,  inelastic  behaviors  lead  to  energy  dis¬ 
sipation.  This  means  that  the  stress-strain  (or  load- 
deformation)  curve  is  not  a  single-valued  function  but 
forms  a  hysteresis  loop.  Energy  is  absorbed  by  the  mate¬ 
rial  system  under  cyclic  load,  and  the  energy  absorbed  is 
proportional  to  the  area  within  the  hysteresis  loop  [46], 
as  in  Figure  2.5. 
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Figure  2,5;  Hysteresis  Loop  of  An  Inelastic  Body 


V' 


Consequently,  another  measurement  of  damping  called 
spe^'fic  damping  capacity  (46]  can  be  obtained  by  compar¬ 
ing  G  energy  dissipated  (or  absorbed)  of  the 

syst„..  in  a  cycle  with  maximum  strain  energy  stored 
(Ug)j(j3j.  in  the  system  i-’  i  -.ng  that  cycle. 

'I’c  “  ^*^0^  cyc^  ^^s^max  2.21 

In  Appendix  A,  viscoelastic  material  is  shown  to  have  such 
a  hysteresis  loop  under  cyclic  loading,  and  the  same 
expression  for  specific  damping  capacity  is  obtained.  The 
difference  is  that  specific  damping  of  viscoelastic 
material  is  also  a  function  of  frequency,  as  reported 
in  studies  [48,  49].  This  is  because  the  storage  and  loss 
moduli  are  functions  of  frequency. 
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For  small  damping,  the  relationships  between  those 
representations  of  damping  are  given  in  references 
[50,  51], 
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CHAPTER  3 

DAMPING  OF  UNIDIRECTIONAL  FIBER  COMPOSITES 

3.1  Introduction 

The  objective  of  this  chapter  is  to  determine  theo¬ 
retically  the  damping  of  unidirectional  fiber  reinforced 
polymer  matrix  composites.  The  major  damping  mechanism  of 
such  composites  is  the  viscoelastic  behavior  of  the  poly¬ 
mer  and  fibers.  The  analysis  is  carried  out  by  first 
applying  the  concepts  of  balance  of  force  and  equal  strain 
energy  on  short-fiber  composite  model  to  determine  the 
longitudinal  modulus  of  short-fiber  composite.  Then  the 
elementary  mechanics  approach  is  used  to  find  the  modulus 
along  the  loading  direction  as  a  function  of  the  mecha¬ 
nical  properties  of  the  fiber  and  matrix  materials.  This 
is  followed  by  applying  the  viscoelastic-elastic 
correspondence  principle  {45,  52]  to  express  the  mechani¬ 
cal  properties  of  the  composite,  fiber,  and  matrix;  then 
after  the  real  and  imaginary  parts  of  complex  modulus  are 
separated,  the  damping  of  the  composite  can  be  obtained. 
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3»2  Damping  Analysis  of  Unidirectional  Fiber  Composites 
■3.2  Short-Fiber  Composite  Model 

The  short-fiber  composite  model  is  composed  of  a 
finite-length  fiber  and  the  polymer  matrix,  as  in 
Figure  3.1(a), 


(a)  (b)  (c> 

Figure  3,1:  Short-Fiber  Composite  Model 


Figure  3.1(c)  is  the  homogeneous  material  equivalent  to 
the  composite  of  Figure  3,1  (a).  Figure  3.1(b)  is  the 
front  middle  longitudinal  section  view  of  Figure  3.1(a), 
where  d  and  s  are  the  diameter  and  the  length  of  fibei; 
respectively;  D  and  L  are  the  diameter  and  the  length  of 
the  composite  model,  respectively;  and  P  is  interpreted  as 
the  distance  between  fiber  tips  along  fiber  direction. 

The  ratio  of  p  to  s  is  defined  as  R  and  is  interpreted  as 
the  degree  of  discontinuity.  During  the  derivation  of 
Young's  modulus  along  the  fiber  direction,  the  short-fiber 
composite  model  is  treated  as  if  it  is  composed  of  two 
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materials  connectedin  series  along  the  fiber  direction. 
One  material  which  is  between  sections  K  and  H'  is  the 
mixture  of  fiber  and  matrix  having  length  s,  while  the 
other  material  is  just  the  pure  matrix  having  length  P. 

As  In  some  other  analytical  work  on  short-fiber 
composites  [23/  31,  32],  the  results  of  Cox's  shear  lag 
stress  analysis  [20]  are  used  in  this  study.  The  expres¬ 
sion  for  elastic  stiffness  of  the  discontinuous  fiber 
composite  is  derived  from  the  average  of  fiber  stress 
based  on  Cox's  fiber  stress  distribution  (in  which  the 
longitudinal  fiber  stress  is  a  function  of  position) . 


nr  =  t-E-  -  cosh[6.(s/2-x)  3  I 
i.  -L  cosh(6s/2)  ' 


where  x,  6,  and  6s/2  are  defined  in  reference  [31],  and 
is  the  strain  of  the  fiber.  In  this  study,  the  square 
packing  array  of  fiber  composites  is  considered;  there¬ 
fore,  6s/2  can  be  written,  according  to  reference  [31],  as 


6s  _  n  s 
2~  "  ^  d 


(  - S - ) 

^f  Ln  ^ 
4Vf 


3.2 


The  average  fiber  stress  is 


o.c  dx 


3.3 
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Substitute  Equation  3.1  into  Equation  3.3 


Uf  =  ef  Ef  [1 


tanh (6s/2)  i 
63/2  ^ 


3.4 


For  the.  composite  between  sections  K  and  H'  in  Figure 
3.1(3),  in  order  to  have  static  equilibrium,  the 
total  longitudinal  force  q  applied  to  this  composite 
must  be 


q  =  =  Of  Af  + 

Therefore, 


'^f  ^f 


3.6 


f  9 

where  and  Vjj,  are  the  fiber  volume  fraction  and  matrix 
volume  fraction  within  sections  K  and  H*  separately. 

It  is  assumed  that  the  composite,  fiber,  and  matrix  (all 
between  sections  H  and  H')  have  the  same  extensional 
strain  e.  The  longitudinal  modulus  of  material  between 
sections  H  and  H’  can  be  obtained  from  Equation  3o6. 


v'  [1  -  tanMg|/2l3 


+  E 


m  ra 


3.7 
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$  f 

If  Vf  and  Vjp  are  expressed  in  terms  of  V^,  and  R,  Equa¬ 

tion  3.7  can  be  rewritten  as 

a  Ef(V^  +  VfR)  [1  -  (Vj^  -  VjR) 

3.8 

Alternatively,  if  the  same  assumption  is  used  as  in  conti¬ 
nuous  fiber  comoosites  is  considered,  the  fiber  stress 
along  longitudinal  direction  is  assumed  to  be  uniform 
everywhere  in  fiber,  and  the  longitudinal  modulus  of 
material  between  sections  H  and  H'  can  be  obtained  by 
using  rule  of  mixtures. 

=  Ej  (V£  *  VjR)  +  E„(V„  -  Vf  R)  3.S 

Equations  3.8  and  3,9  show  that  for  continuous  fiber  com¬ 
posites,  the  longitudinal  Young's  modulus  obtained  by  the 
rule  of  mixtures  is  higher  than  that  obtained  by  Cox's 
analysis.  This  is  because  in  the  rule  of  mixture 
approach,  uniform  longitudinal  fiber  stress  is  assumed, 
while  in  the  Cox's  approach,  uniform  longitudinal  fiber 
stress  exists  only  at  the  locations  far  away  from  the 
fiber  tips,  and  this  longitudinal  fiber  stress  reduces  to 
zero  at  fiber  tips.  Finite  element  stress  analyses  [25, 

26]  show  that  the  reduction  of  long.itudinal  fiber  stress 
around  fiber  tips  does  exist;  and  the  magnitude  of  this 
stress  is  not  zero  but  finite.  So  it  is  hard  to  say  which 
approach  (rule  of  mixtures  or  Cox's  analysis)  is  more 
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nearly  correct.  However,  Table  3.1  shows  the  values  of 
tanh  (Bs/2)/(8s/2)  of  graphite  epoxy  and  Kevlar  epoxy  with 
Vf  being  0.7  or  0.4.  This  table  indicates  that  the  modi¬ 
fication  term  tanh(6s/2) / (Bs/2)  becomes  important  when 
fiber  volume  fraction  and  fiber  aspect  ratio  are  both 
small.  On  the  other  hand,  when  the  fiber  volume  fraction 
is  greater  than  0.4  and  the  fiber  aspect  ratio  is  greater 
than  100,  the  effect  of  tanh(Bs/2) {Bs/2)  could  be 
neglected. 

Table  3.l!  Influence  of  Vf  and  s/d  on  tanh  (Bs/2)  (6s/2) 


Composite 

i 

Vf 

s/d 

tanh  (Bs/2) / (6s/2) 

Graphite-epoxy 

0.4 

5. 

0.725 

Graphite-epoxy 

0.7 

5. 

0.369 

Kevlar-epoxy 

0.4 

5 . 

0 . 610 

Graphite-epoxy 

0.4 

25. 

0.180 

Graphite-epoxy 

0.7 

25. 

0.074 

Kevlar-epoxy 

0.4 

25 . 

0.135 

G  raph i te- epoxy 

0.4 

100. 

0.045 

Wnen  the  same  external  stress  a  is  applied  to  short- 
fiber  composite  model  and  to  its  equivalent  homogeneous 
materials,  the  same  strain  energy  density  U  is  presumed 
for  those  two  materials  (Figure  3.1(a)  and  Figure  3.1(c)) 


3.10 
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and 


2  2 

n-l2_£  +  i  1L-S_  3.11 

a  ■  2  L  2  E.,  L 

U  .  -1-  2^  2.12 

c  2 

where  Ej^  is  the  longitudinal  Young’s  modulus  of  the  homo¬ 
geneous  material  equivalent  to  short-fiber  composite 
model.  After  Equations  3.11  and  3.12  are  substituted  in 
Equation  3.10^  Ej^  can  be  expressed  as 


£  E 

,  _  c  m _ 

Xi  g  ^  g  1 

c  l-»-R  m  1+R 


3.13 


where  E^.  can  be  obtained  by  Equation  3.8  or  Equation  3.9, 
and  R  is  the  ratio  of  P  to  s.  When  R  equals  zero,  Ej^ 
equals  E^.  On  the  other  hand,  if  R  is  a  very  large 
number,  Ej^  will  be  very  close  to  Ejj^.  Since  the  fiber 
aspect  ratio  considered  in  this  study  is  between  25  and 
10,000,  unless  specially  mentioned,  the  modified  Cox's 
analysis  (i.e..  Equations  3.8  and  3.13)  is  used  to 
anaylize  the  longitudinal  Young's  modulus  of  short-fiber 
composite  model. 
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It  should  be  noted  that  the  continuous  fiber  compo¬ 
site  can  be  induced  either  by  letting  the  fiber  aspect 
ratio  be  a  very  large  number  in  the  modified  Cox's  analy 
sis  or  by  letting  R  be  zero  in  the  modified  rule  of 
mixture. 

3.2.2.  Damping  of  Aligned  Short-Fiber  Composites 
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For  a  continuous  aligned  composite,  the  off-axis  modulus 
Ejj  along  the  loading  direction  is  given  in  reference  (381 


cos  0 


4 

sin  0  ^  ^  1 

Eij,  CjL-j 


2  2 
sin  0  cos  0 


3.14 


where  Ej^  and  represent  the  moduli  along  and  transverse 
to  the  fiber  direction  respectively,  is  the  in-plane 

shear  modulus,  and  is  the  major  Poisson's  ratio. 
Equation  3.14  can  be  easily  derived  from  the  elementary 
mechanics  approach  for  off-axis  continuous  fiber 
composites.  For  off-axis  aligned  short-fiber  composites, 
one  can  derive  a  similar  expression  for  Ej^  from  Equation 
3.14  by  replacing  E^,  E^,  G^ij,  and  by  the  correspond¬ 
ing  formula  for  aligned  short-fiber  composites.  The 
longitudinal  modulus  can  be  obtained  from  Equations  3.8 
and  3.13.  The  transverse  modulus  Ej,  in-plane  shear 
modulus  Glt  and  the  major  Poisson's  ratio,  those  material 
constants  are  assumed  to  be  independent  of  length  of 
fiber,  can  be  obtained  by  using  the  Halpin-Tsai  Equation 
[53]  and  the  rule  of  mixtures,  i.e. 


E 

T 


E 

m 


1  +  2n  V 
1  f 


1  -  n^Vf 


3.15 


1  + 


G  »  G 
LT  m 


1  -  n2Vf 


3.16 
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'I.T 


=  V 


fLT 


3. 17 


where 


n 


1 


(E 


fT 


(E 


fT 


/E  )  -  1 
m _ 

^  2 


3.18 


_  fLT  m _ 

2  (GfLT/^m^  +  ^ 


3.19 


Up  to  now,  ail  equations  presented  in  this  chapter 
are  derived  for  elastic  material.  When  a  viscoelastic 
material  is  considered,  the  elastic-viscoelastic 
correspondence  principle  (52]  can  be  used  to  obtain  the 
corresponding  relationships  of  viscoelastic  material.  For 
viscoelastic  material  the  basic  material  properties  are 


redefined  as 
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i  ^fL 

E|t  -  +  i 

<^fLT  “  ^kx"^  ^  ^kT 


E* 


EA 

m 


+  i  E 


m 


Gm 


m 


+  i 


G“ 

m 


m 


=  V ' 


m 


+  i  V 


m 


'^fLT 


'fLT 


3.20 


E*  s  E*  +  i  E"  3.21 

^  X  X 

Where  the  prime  quantities  indicate  the  storage  moduli  or 
storage  Poisson's,  the  double  prime  quantities  indicate 
the  loss  moduli  or  loss  Poisson's  ratio,  and  the  i  is 
defined  as  FI.  In  this  research,  bulk  modulus  of  epoxy 
matrix  Kj^  is  assumed  to  be  real  and  Independent  of  fre¬ 
quency  (151.  For  Isotropic  epoxy  matrix. 


nt 


m 


rrr^ 


m 


3.22 
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While  the  viscoelastic  behavior  of  epoxy  is  considered, 


\  = 


TTT^ 


K  +  E" 
m  m 

'  "  i 
m 


m 
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The  complex  form  of  Vjjj  can  be  obtained  by  Equation  3.23 


V '  +  i  v"  = 
m  m 


,  E •  +  i  E" 

2  ■  — nr - ^ 

m 


3.24 


4>^Cvsji»  <3  C4  V*4>W44  \/Jm  Wii^  i.coi«>  cai&SA  wsuli  vi^  a.iA  wka*  ^  ^4ifc.Li»o  wJ. 

hand  side  of  Equation  3.24,  one  will  have 


v'  +  i  v’ 
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m 


m 


+  i  (v;  -  |)  ] 
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Similarily,  the  complex  form  of  shear  modulus,  0'^^^+  G"^  of 
viscoelastic  matrix  can  be  expressed  as  functions  of  E^, 


Em»  i-®-' 


G' 

m 


K  K 

.  1  t-H  „ _ J.  ^ 

m  2  d+v' )  I’d+V ) ' 

m  m 
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The  complex  form  of  6s/2  of  Equation  3.3,  as  shown  in 
Appendix  B,  is 


G**  E  * 

6*s  .  i  B*s  _  6s  .  6s  ,_m  _  3.27 

"*2  ^4  'g*  E*  ^ 

m  f  L 


The  reason  that  the  imaginary  part  of  fiber  Poisson's 

ratio  is  set  as  zero  in  the  last  equation  of  Equation  3.20 

is  because  first,  most  fibers  are  known  to  be  anisotropic 

materials;  therefore.  Equation  3.22  is  not  true  for  most 

fibers.  Secondly,  the  corresponding  term  of  for 

most  fibers  (i.e.,  /E^) ,  except  Kevlar  fiber,  is  much 

less  than  E"/E'  .  Due  to  the  lack  of  the  available  data 
m  m 

of  transverse  damping  and  shear  damping,  those  two  dampings 

are  assumed  to  be  equal  to  E"  /F'  . 

£L  £L 

It  should  b,  noted  that  n  and  n,  are  treated  as 

ro  f 

as  material  properties,  and  they  are  defined  as 
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E" 


-  ^ 

m  ”  E  * 
m 


3.28 


n 


"■fL 
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By  using  Equations  25  and  3.28-3.29,  the  right-hand 
side  of  Equations  3.20,  3.21,  and  3.27  can  be  rewritten  as 
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i  Ecm  “  (1  +  i  n.  ) 

fT  fT  fT  f 


+  i  G"  *  G'  ^(1  +  i  n,  ) 
fLT  fLT  fLT  f 


E*  +  1  E" 
m  m 


=  E’  (1  +  i  n  ) 
ra  m 


3.30 


G'  +  i  G"  =  G'  (1  +  i  n_  ) 
ni  m  m  Gm 


v'  +  i  v" 
m  m 


v'  +  i  n  (v* 
m  mm 


1\  * 
i)  =  V 


V '  -  V 

fLT  fLT 


E'  +  i  E" 

X  X 


^  i  lls  =  is  {}  ->-  i  i  (f)  -n)l  = 

2  “  2  2  “  2  Gm  f  '  2 


where  ri  is  defined  as 
Gm 


'^Gm 


3.31 


From  Equation  3.31,  it  is  observed  that  rigm  higher  than 
After  using  Equations  3.20,  3.27,  and  3.30,  one 
can  rewrite  Equations  3.13  and  3.15-3.17  for  viscoelastic 
material,  as  follows 
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(e:  +  i  e:)  (e^  +  1  e;) 
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<E’  +  i  E-)  +  (E^  +  e;;) 


3.32 


(El  +  i  E") 
m  m 


1  +  2nJ 
i  -  Vf” 
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;*  =  {G'  +  i  G-) 

LT  m  m 


1  - 
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V*  =  (v'  +  V  V  )  +  i  n  (v'  -  ■i)V 

LT  fLT  f  mm  m  m  2  m 


3.35 


where 


E'  +  i  E" 
c  c 


[1  -  tanh(^6s/2)  j  +  +  iE"jjj )  V^R) 

6  s/2  "*  3.36 


tanh{B*s/2)  =  tanh{6s/2) 


n  n, 

+  i  M 


cosh  (Bs/2) 


3.37 
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n?  = 


<^;lt  ^  ^  ^  i 


n\ 
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After  substitute  e£,  e£,  G^t  /  and  v£,j,  from  Equations 
3.32-3.35  for  E^/  E.J.,  G,j,,  and  respectively,  and 

Ej^+iEjj  for  Ejj  into  Equation  3.14,  one  obtains 
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IPTIE^ 
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cos^  e  .  sin*^  0  ^  ,  1 

— X -  +  - 3f -  ' — sr- 
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i£)cos^0  sin^0 


3.40 


Damping  of  the  aligned  short-fiber  composite  along  x 
direction,  rij{»  is  then  determined  by 


K 

■"E^ 

X 


3.41 


Equations  3.40  and  3.41  show  that  material  damping  and 
stiffness  of  aligned  short-fiber  composite  are  functions 
of  material  properties  of  fiber 

VfTT/  and  n^)  and  matrix  (i.e.,  E',  V,  and  n  ) ,  fiber 
rjj-i-  f  m  m  m 

aspect  ratio  (s/d),  fiber  volume  fraction  (Vj),  degree  of 
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discontinuity  (R) ,  loading  direction  (6),  and  packing 
geometry  of  fiber.  If  the  four  different  kinds  of  pre¬ 
packed  tapes  (glass-epoxy,  Kevlar-epoxy,  graphite-epoxy, 
and  boron-epoxy)  are  used  to  make  the  fiber  composite,  the 
design  variables  utilized  to  analyze  the  stiffness  and 
material  damping  are  s/d,  0,  R,  ^fT*  ^fLT^ 

''fLT  '  ''f  '  Nn'  'I 


(E^,  s/d,  e,  R,  V^)  3.42 

=  f2  (Ef/  E^/  s/d,  e,  R,  Vj)  3.43 

A  similar  approach  can  be  applied  to  determine  damping 
along  y  direction,  r\y,  and  stiffness  along  y  direction. 


CHAPTER  4 

DAMPING  OF  RANDOMLY  ORIENTED  SHORT-FIBER  COMPOSITES 

4.1  Introduction 

The  objective  of  this  chapter  is  to  determine 
analytically  the  material  damping  of  in-plane  randomly 
oriented  short-fiber  composites-  The  analysis  is 
carried  out  by  using  the  extension  of  the  short-fiber  com¬ 
posite  model  and  part  of  the  results  obtained  in 
Chapter  3.  An  averaging  procedure  is  first  applied  to 
the  six  off-axis  reduced  stiffnesses  (i,  j  =  1#  2,  6) 

with  respect  to  the  angle  e  between  the  fiber  orienta¬ 
tion  and  the  applied  load.  The  results  of  integration 
show  that  in-plane  randomly  oriented  fiber  composites 
behave  like  a  planar  isotropic  material.  By  using  the 
properties  of  isotropic  materials.  Young's  and  shear 
moduli  can  be  obtained  as  functions  of  the  reduced 
stiffnesses  Qj^j  (1,  j  =  1,  2,  6).  After  the  application 
of  the  elastic-viscoelastic  correspondence  principle  and 
separation  of  the  real  and  imaginary  parts  of  the  complex 
Young's  and  shear  moduli,  material  damping  is  obtained. 
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4.2  Dampinq  Analysis  of  in-Plane  Randomly  Oriented 
Short-Fiber  Composites 

For  in-plane  randomly  oriented  short-fiber  compo¬ 
sites,  no  difference  caused  by  different  direction  paral¬ 
lel  to  the  planes  on  which  fibers  are  laid.  The  averaging 
procedure  i  "jne  of  the  approaches  which  will  lead  to  the 
Isotrop.  Therefore,  an  averaging  by  integrating  the  six 
moduli  of  off-axis  short-fiber  composites  with  respect  to 
0  from  0*0  to  0=Tt  should  be  used.  However,  from  Equation 
3.14  for  Ejj  and  similar  Equations  for  Ey,  ^xy'  ''xy'  "^x 
my  i38],  one  finds  that  it  is  not  convenient  to  integrate 
and  obtain  the  average  in  closed  form  in  terms  of 
Ej^,  E,j.,  and  Instead  of  integrating  the  six 

engineering  moduli,  one  can  integrate  the  six  components 
of  the  off-axis  reduced  stiffness  of  the  plane  stress 
case  (i,  j  =  1,  2,  6)  and  obtain  the  average 
Q j #  i.e. 

Qij  =  ^  [  Qij  i/  j  =  1/  2,  6  4.1 

^  o 

The  expression  for  Q_j^j  (i,  j  =  1,  2,  6)  as  a  function  of  0 
can  be  found  in  reference  138].  After  integrating  with 
respect  to  G  from  0=o  to  G=ii  and  then  dividing  each  of  the 
result  by  ti,  one  obtains 
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Qll  *  ^22  *  i  <Qll  ^22^  +  Qi2^ 

^66  *  g  ^^11  *  ^22^  “  ?  ®12  ■*■  2  ^66 
Qi2  =  g  ^Qll  ■*■  ^22^  J  ^12  ”  i  ^66 
^16  =  Q66  =  0 


It  is  easy  to  show  from  Equations  4. 2-4. 4  that  the  follow¬ 
ing  relation  exists 


Qi2  ■*■  2^66  “  ^11 

Therefore,  after  integration,  there  are  only  two  inde- 
pendent  material  constants,  namely,  Qj.  and  G^. 


Qr  =  Qll 

Gj  =  Bos 

This  implies,  as  expected,  that  in-plane  randomly  oriented 
short-fiber  composites  behave  as  planar  isotropic  mate¬ 
rials  wi  h  two  independent  material  constants.  The 
subscript  r  represents  randomly  oriented  short-fiber 
composites. 


\.*u  ^  »VV  ^  ^  \  /VV  ^  1  ^  \  V  ^  V  *>■  \  ^  T»  ^  ^  \.rm  ^  L  »  -  ^  -w.  > 
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For  isotropic  materials,  the  following  relations 

exist 


mm 

Q 


r 


4.9 


G  =  _ _  4.10 

2(1  +  V) 

where  *v  is  defined  by  Qi2/Qli*  Elimination  of  v  from 
Equations  4.9  and  4.10  yields  the  expression  for 
E^,  the  Young's  modulus  of  a  randomly  oriented  short- 
fiber  composite  as  a  function  of  G^.  and  Q^,. 


Zp  =  4  Gj.  (1 


4.11 


Substitution  of  Equations  4.2,  4.3,  4.7,  and  4.8  in  Equa 
tion  4.11  yields  E^^  as  a  function  of  the  four  reduced 
stiffness  Qii,  Q22'  ^12'  ^66* 


^r  =  (Qll  +  Q22)  -  Q12  2Q66^ 


II  - 


Qll  ^  Q22  ’  20^2  ^  ^066 
"3  (Q^I  +  Q22)  +  2  (Qi2  +  2Qgg) 


4.12 


Similarly, 
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"  I  ^22^  “  1^12  1^66 

Since  Qj^x'  ^22'  ^12'  ^66  directly  related  to  the 

four  basic  engineering  constants  Ej^,  E,j,,  ^LT'  '^LT 

(defined  in  Equations  3.13,  3.15,  3.16,  and  3.17),  accord- 
ing  to  Equations  4.12  and  4.13,  E^  and  Gj.  can  be  expressed 
as  functions  of  E^,  E,j,,  0^^'  '^LT* 

Next,  as  in  Section  3.2.2  for  aligned  short-fiber  com¬ 
posites,  according  to  the  elastic-viscoelastic  correspond¬ 
ence  principle,  one  may  replace  E^  by  E*  =  +  i  2^,  G^. 

by  G*  =  G'^  +  iG^*  ,  E^  by  e£  ,  E,j,  by  ,  C^ip  by  g£,ji  , 

v^T  ^LT  '  ^LT  '  ''lt  '  defined  in 

Equations  3.32-3.35.  After  separation  of  the  real  and 
imaginary  parts,  the  material  damping  constants  tIj.  and  Hqj- 
of  in-plane  randomly  oriented  short-fiber  composites  can 
be  obtained. 


E" 


4.15 


CHAPTER  5 

DAMPING  OF  LAMINATED  FIBER  COMPOSITES — 

LAMINATED  PLATE  THEORY  APPROACH 

5 -I  Introduction 

In  this  study,  laminated  plate  theory  and  an  energy 
approach  ate  used  to  analyze  the  material  damping  and 
stiffness  of  symmetrically  laminated  fiber  composites.  In 
this  chapter,  we  will  discuss  all  analytical  work  of 
laminated  plate  theory  approach,  while  in  Chapter  6  the 
energy  approach  will  be  presented. 

According  to  laminated  plate  theory,  the  constitutive 
equations  (Equation  5.1)  have  already  been  given  in  refer¬ 
ences  (38/  50],  in  terras  of  [A],  [B],  and  [D]  (l.e.  [A]*, 

IB]*,  and  [D] *)  matrices.  Material  damping  of  laminated 
composites  can  then  be  derived  from  the  expression  of  [A], 
[B] ,  and  [D] . 

5.2  Damping  Analysis  of  Laminated  Fiber  Composites 
Through  Lcuninated  Plate  Theory  Approach 

For  a  Icuninated  fiber  composite  plate,  as  in  Figure 
5.1,  the  constitutive  equations  are  given  in  references  [38, 
50]  as  shown  in  Equation  5.1. 
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In  Equation  5.1,  e^,  e^,  and  Yxy  middle  plane  strains, 
k„,  k  and  k„„  are  plate  curvatures  defined  in  [38,  50], 
and  Aj^j,  Bj^j,  and  D^j  (i,  j  =*  1 ,  2,  6)  are  the  equivalent 
reduced  in-plane  stiffness,  coupling  stiffness  and 
reduced  flexure  stiffness,  respectively.  They  are  expressed 
in  terms  of  Q^j  and  total  thickness  of  the  plate,  h,  as 
fol lows ; 


‘ij 


r  h/2  _ 

Qij  dz 

-h/2 


n 


k=l 


5.2a 


Bij 


h/2 


-h/2 


Z  dz  = 


n 


icfl  5.2b 


h/2 


Dij 


dz  = 


n 


-h/2 


5.2c 
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where  Qj^j  can  be  expressed  in  vector  form  [50] 


and 

(i  =  1/  2  .  .  5)  are  defined  in  [50]  as 

U2  =  I  (4Qij^  -  4Q22) 

*^3  ®  ■*■  ^22  "  "  ^^66^ 

^4  “  5  ^Qll  '*'  ^22  ^*212  "  ^^66^ 

i’5  =  g  (Qii  +  Q22  ~  2^12  ■*■ 
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Qij  (i/  j  *  1/  6)  are  known  to  be,  by  references  138, 

501  , 


Q  =  - - 

11  '^TL 

— 

22  ^“'^LT  '^TL 

^  5.5 

_  ^LT 

12  I-'^LT  ^TL 

^66"  ‘^LT 

and 

-  ^LT  5.6 

TL  El 

where  El»  ^LT'  ^LT  defined  in  Equations  3.13, 

3.15-3.17.  Interesting  relations  should  be  noted  that 
Qj^j  (1,  j  =  1,  2,  6)  of  Equation  5.3  will  be  equal  to 
of  Equations  4. 2-4. 5,  providing  U2  and  U3  are  set  as  zero. 
Equation  5.1  can  be  rewritten  in  matrix  form  as 


• 

f 

•  ’ 

N 

A 

B 

e' 

< 

)■  = 

> 

M 

B 

k 

D 

4 

k 

.  > 

lO 
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For  symmetric  leunlnates,  the  coupling  stiffness  matrix  B 
Is  a  zero  matrix,  and  Equations  5.1  and  5.7  can  be 
uncoupled  as 


or 


and 


or 


^x 

^11 

Ai2 

^16 

N 

>  = 

Aj^2 

^22 

^2b 

i 

^  ■ 

^xy 

^16 

^26 

^66 

^xy 

{n} 

= 

[A] 

{eM 

'°11 

Di2 

Die' 

►  = 

^12 

^22 

D26 

< 

ky  > 

Mxy 

.Die 

^26 

Dee 

^xy 

{m}  -  [D]  {k} 


5.8a 


5.8b 


5.9a 


5.9b 


From  Equation  5.8a  and  the  definition  of  one  can 

express  in-plane  moduli  (i,  j  =  1,  2,  6)  as 


E 

Ij 


5.10 
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where  h  Is  the  total  thickness  o£  laminated  composite  and 
are  the  elements  of  the  Inverse  of  matariic  [A^j],  and 
A^j  are  determined  by  the  following  system  of  equations. 


-1  -1  -1 

All  ^12  ^16 

-1  -1  -1 

*  ' 

*11  *12  *16 

'  ' 

10  0 

^12  *22  *26 

-1  -1  “1 

*12  *22  *26 

0  10 

.*16  *26  *66. 

.*16  *26  *66^ 

0  0  1 

5 

According  to  elastlc-vlscoelastic  correspondence 
principle,  when  the  viscoelastic  behaviors  of  fiber 
composites  are  considered  the  elastic  material  constants 
S,p,  G^j,,  and  (defined  in  Equations  3.13,  3.15-3.17) 

should  be  replaced  by  the  corresponding  complex  moduli, 
respectively.  Consequently,  the  complex  form  of  Equations 
5.8  and  5.9  can  be  written  as 


*  1 

Nx 

’ 

*11 

* 

*12 

♦  ^ 

^16 

*  « 

‘'X 

* 

*12 

*22 

*26 

< 

.  "xy. 

*16 

* 

*26 

*66. 

Yxy 

5.12a 


or 
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{n  } 

9 

[A 

]  {eM 

(  *  '' 

★ 

* 

* 

Mx 

^11 

J^12 

*^16 

* 

* 

* 

* 

\  = 

^12 

*^22 

^26 

* 

★ 

* 

* 

^xy 

«16 

°26 

^66 

t  s 


xy 


5.12b 


5.13a 


or 


{m  }  =  [D  ]  {k} 


5.13b 


The  corresponding  complex  form  of  in-plane  complex  moduli 
E*j  can  be  expressed  as 


1  »  " 


5.14 


(^ij) 


*  -1 

where  (A.^)  is  the  element  in  the  ith  row  and  jth  column 
of  the  inverse  of  matrix  [A^j]. 

The  in-plane  material  damping  and  fjLi^^are 

material  damping  are  defined  as  follows: 


n 

I  ij 


iii 


E'  . 
13 


i,  j  =  1,  2,  6 


5.15 
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F^j 


I,  j  -  1,  2,  6 


5.16 


For  the  same  kind  of  fiber  composite,  the  result 
obtained  by  Equation  5.15  for  sixteen  plies  of  unidirec¬ 
tional  laminated  composite  is  same  as  that  obtained  by 
Equation  3.41  for  laminar  composite  with  the  same 
off-axis  angle  ©.  This  indicatv,s  that  the  approach 
presented  in  this  chapter  is  correct,  although  it  may  not 
be  convenient  because  the  inverse  of  a  complex  matrix  is 
involved. 


CHAPTER  6 

DAMPING  OF  LAMINATED  FIBER  COMPOSITES—ENERGY  APPROACH 

6.1  Introduction 

The  drawback  of  damping  analysis  using  laminated 
plate  theory  approach  is  that  it  does  not  include  the 
effect  cf  interlaminar  stresses  (the  stresses  at  the 
interfaces  of  laminated  composites).  It  has  been  shown  in 
study  [54]  that  even  when  an  in-plane  uniform  tension  load 
is  applied  to  the  laminated  composites,  there  do  exist 
appreciable  interlaminar  stresses  around  the  free  edges. 

In  this  chapter,  an  energy  approach  in  conjunction 
with  a  three-dimensional  finite-element  method  [55,  56]  is 
used  to  analyze  the  material  damping  under  certain  loading 
and  boundary  conditions.  It  is  believed  that  this  model 
represents  a  more  realistic  approach  by  including  the 
energy  dissipated  at  the  interfaces.  By  using  this 
approach,  we  improve  the  approach  from  a  two-dimensional, 
classical,  laminated  plate  theory  to  a  three-dimensional 
elasticity  theory. 
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The  damping  of  laminated  materials  in  the  first  mode 


vibration  is  determined  as 


n 


^  eye 

q=l _ _ _ 

n 

Z  2Tr  Ug 
q=l 


6.1 


where  n  is  the  total  number  of  the  plies,  (q^D^r-yc 
energy  dissipated  in  the  qth  layer  during  a  cycle,  and 
is  the  maximum  strain  energy  stored  in  the  qth  layer. 
Detailed  expressions  of  (qUo)cyc  q^s  given  in 

Appendix  A  in  which  the  energy  expressions  for  a  visco¬ 
elastic  material  given  in  reference  157]  are  used 

The  analytical  expression  of  the  maximum  strain 
energy  Ug  for  an  elastic  body  is 


(j,  k  =  1,  2  .  .  .  6) 


where  Cj  are  the  maximum  strains,  and  Cjj^  are  the  moduli 
of  the  elastic  body. 

The  expressions  of  Cjj^  (j,  k  =  1,  2  ...  6)  for 
orthotropic  elastic  material  are  given  by  Jones  [58]; 
here  k  and  j  axes  are  material  principal  axes.  In  this 
study,  each  layer  of  fiber  composite  is  considered  to  be 


56 


transversely  isotropic  material.  Therefore,  as  presented 
in  Appendix  B.l,  the  moduli  of  mth  layer  fiber  composite 
could  be  defined  as  functions 


^^jk^m  -  f  ,  E-r  ,  Gl,j,  , 


where  ,  E,j,  ,  ,  and  Vj^^,  are  given  in  Section  3.1  for 

short-fiber  composites,  and  Vq,.p,  is  assumed  equal  to  v^. 

An  analytical  expression  for  is  available  in  the 

literature  [38],  but  the  assumption  that  =  Vjj^,  is 

believed  to  be  accurate  enough  in  our  analysis.  According 
to  the  elastic-viscoelastic  correspondence  principle,  when 

the  viscoelastic  behaviors  of  the  body  are  considered,  Ej^, 

_*  *  *  * 

Eiji,  should  be  replaced  by  E^,  hj,, 

* 

and  respectively.  For  example. 


''tl  ''lt 


changes  into 


M  •  m"  P  <  P" 

(  TL  +  i  TL)  (  L  +  i  L) 

Tv"* - — +~rv  ""T 

^  TL  '  LT  ^  LT' 


^12  "  ^12  ^  ^"l2 


t ..  1.  V.  Vn  I  \  W. V JV'\.S.'\.'v"Jv\v'.'V\'v\  -s.'  -."i  a.V  ".r 
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After  separating  the  real  and  imaginary  parts  and  using 
matrix  rotation  (if  the  material  principal  axes  do  not 
coincide  with  the  global  axes),  the  complex  moduli  C* 
with  respect  to  the  global  axes  (i.e.,  x,  y,  and  z  axes] 
of  the  kth  layer  of  fiber  composites  are  then  obtained  in 
the  form 

=  tC'  ]j^  +  i  6.7 

where  C  are  storage  moduli,  and  C"  are  loss  moduli.  All 
three  matrices  are  symmetric  matrices. 

The  storage  energy  and  the  dissipated  energy 
cinrinrr  a  f  T!  ^  In  the  oth  1  aver  of  laminated 

composites  can  be  expressed  as 

Cjk  ^k  6.8 


^jk  ^k 


q^s 


1 

2 


and 


<qUD)cyc  “  ^ 


correspondingly.  Consequently,  the  material  damping  of 
a  n  layers  laminated  composite  in  the  first  mode  vibration 


58 


Z  I  {e}’^  (C-]  {e}  dv 
q=l  ■'  q 

■■  ■  '  - ,  — — - 

Z  [  {e}'^  {C']  (e}  dv 

q»l  Jvq 


6.10 


where  strain  field  {e}  depends  on  the  loading  and  boundary 
conditions,  and  the  strain  field  would  be  determined  in 
this  study  through  a  three-dimensional  finite  element 
method. 

From  Equation  6.10,  one  may  thus  arrive  that  even  a 
highly  dissipative  modulus  cannot  contribute  significantly 
to  the  total  loss  factor  n#  if  it’s  associated  strain  (or 
stress)  does  not  participate  considerably  in  the  total 
stored  energy. 

The  procedures  taken  in  energy  approach  are  briefly 
described  as  follows: 

Step  It 

As  in  the  laminated  theory  approach,  the  elastic 
solution  is  first  sought  in  energy  approach.  The 
equation  of  motion  of  an  elastic  body  in  static  case  can 
be  derived  from  the  principle  of  stationary  potential 
energy  given  in  reference  [59] . 

6(U-W^)  =  0  6.11 


where  U  and  are  the  strain  energy  and  the  work  done  by 
external  forces,  respectively. 
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Step  2; 

After  the  assumed  displacement  function  (as  a 
function  of  nodal  displacement  q)  is  substituted  into 
Equation  6.11,  the  equations  of  motion  can  be  expressed  as 
a  system  of  equations  (details  are  given  in  Appendix  B.2). 

[k]  fi;  =  {?}  6.12 

where  [k]  is  the  stiffness  matrix,  and  {f}  is  the  nodal 
forces . 

Step  3 : 

After  the  displacement  field  is  determined  by  sub¬ 
stituting  the  solution  of  Equation  6.12  into  the  assumed 
displacement  function,  the  strain  field  (e}  of  the  elastic 
body  could  be  obtained  through  displacement-strain 
relations . 

Step  4; 

Once  the  strain  field  of  an  elastic  body  is  known,  the 
material  damping  of  a  viscoelastic  body  under  the  same 
loading  and  boundary  conditions  can  be  determined  by 
Equation  6.10. 

In  general,  we  can  apply  either  uniform  normal  and 
uniform  shear  forces  (N^^,  Ny,  or  N^y)  at  edges  x  (or  y)  = 
constant  or  apply  uniform  normal  and  uniform  twisting 
moments  at  x  (or  y)  =  constant.  In  the  former  cases,  we 
obtain  in-plane  normal  and  shear  damping  .t'^66^' 

and  in  the  latter  cases,  we  obtain  flexural  normal  and 
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twisting  damping  (pOxi»  f*^66^*  Only  the  case  of  uniform 
Njj  applied  at  x  =  constant  will  be  presented  for 
different  symmetric  cross-ply  laminates  in  Chapter  9. 
Other  loading  conditions  can  also  be  accommodated  in  this 
approach. 


CHAPTER  7 

EXPERIMENTAL  MEASUREMENT  OF  DAMPING 

7.1  Introduction 

Free  vibration  decay  {91,  band-width  method  [9], 
resonant-dwell  method  (601,  forced-vibration  techniques 
[61],  and  impulse  techniques  [62,  63,  64]  are  the  very 
popular  experimental  techniques  used  to  measure  the 
material  damping.  All  of  these  tests  are  subject  to  the 
air  drag  [65,  66],  if  the  tests  are  not  conducted  under  a 
vacuum  condition.  In  this  study,  an  improved  impulse 
technique  approach  [64]  is  utilized  to  measure  the 
material  damping. 

The  impulse  technique  consists  of  the  application  of 
a  force  pulse  at  a  point  on  the  test  structure  and  the 
measurement  of  the  response  at  another  point.  The  input 
force  and  response  signal  are  digitally  processed  by  the 
analyzer  to  form  the  frequency  response  function  or 
transfer  function.  Damping  and  natural  frequencies 
can  then  be  extracted  from  the  output  of  the  analyzer. 

All  measures  are  performed  using  a  digital  signal  process¬ 
ing  technique.  It  involves  filtering  and  sampling  the 
input  wave  forms.  The  sampling  process  converts  a 
voltage  (at  a  certain  point  in  time)  into  a  numerical 
representation. 
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These  numbers  are  then  processed  digitally  to 
produce  tiie  various  calculations  performed  by  the 
analyzer.  The  primary  output  from  the  Fourier  analyzer 
is  the  frequency  response  function,  which  is  a  measure  of 
system’s  characteristics.  The  analyzer  also  calculates 
the  coherence  function,  which  ranges  from  0.0  to  1.0.  The 
larger  the  potential  measurement  noise  or  the  system  non¬ 
linearity,  the  lower  the  coherence  function  will  be  [67], 

A  coherence  value  of  unity  indicates  that  the  output  is 
completely  related  to  the  input.  Thus,  the  coherence 
function  used  here  is  a  measure  of  the  "quality"  of  the 
data. 


7 . 2  Apparatus 

The  composite  plates  or  beams  are  fixed  by  two  alumi¬ 
num  blocks  at  one  end.  A  non-contact  probe  {KD-23X0-3U, 
Kaman  Science  Corporation,  Colorado  Springs,  Colorado), 
known  as  a  motion  transducer,  is  located  about  1.5  mm 
below  the  tip  of  each  specimen.  This  motion  transducer 
operates  under  the  principle  of  eddy  current.  An  aluminum 
foil  target  of  diameter  20  mm  is  cemented  underneath  the 
tip  of  the  specimen.  A  force  transducer  is  mounted  to  the 
head  of  the  impulse  hammer  (Model  K291A,  Piezotronlcs, 

Inc.,  New  York),  to  measure  the  force  input  to  the  speci¬ 
men.  The  other  tip  of  the  hammer  is  connected  to  the 
fixed  spring  so  that  the  magnitude,  location,  and  the 
dwelling  time  of  the  Impact  load  can  be  controlled.  To 
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improve  the  coherence,  an  octave  filter  (4302  dual  24db, 
Ithaco,  Inc.,  New  York),  which  is  connected  to  the  force 
transducer  at  the  hammer  tip  on  one  end  and  to  the  Fast 
Fourier  Analyzer  (FTT,  Model  5420,  Hewlett  Packard)  on  the 
other  end,  is  used  to  magnify  the  input  signal  100  times. 
The  impact  point  is  near  the  stiffer  portion  of  the  speci¬ 
men.  The  excitation  and  response  signals  are  fed  into  the 
FFT  analyzer,  which  displays  the  frequency  response  func¬ 
tion  and  coherence  function.  Each  frequency  response  and 
each  coherence  frunction  are  based  on  a  statistical  analy¬ 
sis  of  an  ensemble  of  six  tests.  A  schematic  drawing  of 
the  experimental  set-up  is  shown  in  Figure  7.1. 

The  frequency  response  function  is  a  complex  valued 
function.  A  typical  experimental  display  of  the  real  and 
imaginary  parts  of  the  frequency  response  function  for  a 
graphite-epoxy  composite,  according  to  reference  [64],  is 
shown  in  Figure  7,2  and  an  enlarged  schematic  drawing  of 
real  part  of  frequency  response  function  is  shown  in 
Figure  7,3.  As  prescribed  in  references  (62,  63,  64],  the 
p?ak  of  the  imaginary  part  determines  the  resonant  fre¬ 
quency,  and  then  from  the  corresponding  real  part  (see 
Figure  7.3)  the  material  damping  (loss  factor,  n)  can  be 
calculated  by 


(ft/fO  V L 


n 


7.1 
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In  using  the  impulse  technique  to  measure  the  mate¬ 
rial  damping,  following  precautions,  as  reported  in  study 
(64],  must  be  taken. 

1.  In  order  to  improve  the  accuracy  of  experimental  data 
of  damping  associated  with  a  particular  mode  of 
vibration,  the  values  of  natural  frequency  within  the 
frequency  range  of  interest  (say  0  to  1,600  Hz)  are 
first  approximately  determined.  Then  a  zooming 
technique  to  increase  the  resolution  of  the  response 
in  the  neighborhood  of  this  particular  frequency  is 
used. 

2.  It  is  necessary  to  avoid  measurement  of  response  near 
a  nodal  point  for  the  modes  to  be  tested.  Such  meas¬ 
urements  would  consist  primarily  of  noise,  since  the 
actural  response  is  very  small  near  nodal  points. 

3.  The  amplitude  of  the  vibration  must  be  kept  below  the 
thickness  of  the  specimen  to  ensure  that  air  damping 
(air  drag)  is  negligible,  since  the  air  damping  is 
linearly  proportional  to  the  amplitude-to-thickness 
ratio  of  the  beam  specimen  [66]. 

4.  It  is  important  to  optimize  the  Analog  to  Digital 
Converter  (ADC)  range  setting  on  an  FFT  analyzer 
before  making  a  measurement,  since  an  optimized  ADC 
range  set  will  increase  resolution  in  the  digitizing 
process  135). 

Recently,  Suarez  and  Gibson  (35,  36],  did  some 
experiments  on  material  damping  of  short-fiber  composite 
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through  impulse  hammer  technique.  Some  of  those  experi¬ 
mental  results  for  unidirectional  short-fiber  composites 
presented  in  this  chapter  are  compared  with  the  analytical 
results. 


Lgure  7,1:  Schematic  Drawing  of  the  Experimental 
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Figure 


7.3;  Enlarged  Schematic  Drawing  of  Real  Part  of 
Frequency  Response  Function 


CHAPTER  8 

OPTIMIZATION  OF  DAMPING  AND  SPECIFIC  STIFFNESS 
FOR  FIBER  COMPOSITES 


8.1  Introduction 

Plunkett  and  Lee  [19]  pointed  out  that  maximum  damp¬ 
ing  can  be  obtained  through  properly  choosing  the  length 
and  spacing  of  the  constraining  layers.  Those  layers  are 
used  to  constrain  the  viscoelastic  layer  (or  coating)  on 
the  structure  surface.  In  the  study  presented  here,  the 


optimisation  of  damping  is  based  on  the  composite  material 


itself  (i.e.,  choosing  the  proper  fiber  length,  fiber 
direction,  the  longitudinal  distance  between  fiber  tips, 
etc.).  It  is  known  that  high  damping  can  reduce  the 
displacement  at  and  near  the  resonant  frequencies,  while 
high  stiffness  can  further  reduce  the  displacement  at 
other  frequencies.  The  best  way  to  optimize  both  damping 
and  stiffness  is  to  keep  them  as  high  as  possible.  Unfor¬ 
tunately,  there  exists  a  general  trend  that  higher  damping 
is  mostly  coupled  with  low  stiffness  and  vice  versa. 
Therefore,  the  idea  of  the  optimization  in  this  study  is 
to  try  to  increase  damping  without  sacrificing  the 
stiffness  too  much. 

From  the  formulations  presented  in  previous  chap¬ 
ters,  one  can  see  the  complication  contained  in  damping 
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analysis.  It  will  be  a  simplier  approach  for  the  optimi¬ 
zation  analysis  if  no  derivative  is  needed.  Consequently, 
the  so-called  Sequential  Simplex  Method  [68,  69]  has  been 
selected  to  analyze  the  optimization. 

8.2  Brief  Introduction  to  Sequential  Simplex  Method 
The  sequential  Simplex  Method  [68,  691  takes  a  regu¬ 
lar  geometric  figure  (known  as  a  simplex)  as  a  base. 

Thus,  in  two  dimensions  (i.'=^.,  two-design  variables),  one 
should  choose  an  equilateral  triangle;  and  in  three  dimen¬ 
sions  (i.e.,  three-design  variables),  one  should  choose  a 
tetrahedron. 

Qbser’.’’ations  (experiments)  are  located  so  that  the 
objective  function  is  evaluated  at  the  points  formed  by 
vertices  of  tht  geometric  figure.  One  vertex  is  then 
re  ted  as  being  inferior  in  value  to  the  others.  The 
r^l  direction  of  search  may  then  be  taken  iu  a  direj- 
.  ay  from  this  worst  point,  the  direction  being 
cl  so  that  the  movement  passes  through  the  center  of 
gravity  of  the  remaining  points.  A  new  point  is  then 
selected  along  this  direction  so  as  to  preserve  the  geome¬ 
tric  shape  ef  the  figure,  and  the  function  is  evaluated 
anew  at  this  point.  The  method  proceeds  with  this  process 
of  vertex  rejection  and  regeneration  until  the  figure 
straddles  the  optimum.  When  no  subsequent  moves  would 
lead  to  further  improvement,  the  last  few  geometric 
figures  are  essenrially  repeated.  There  are  three  rules 
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that  govern  the  whole  procedure  of  this  method-  These 
rules  are  explained  in  the  two-dimensional  case  as 
fol lows  T 

1.  Take  the  point  to  be  rejected  where  the  worst  value 
of  the  objective  function  is  obtained,  and  replace  it 
by  its  reflection  in  the  opposite  side  of  the  tri¬ 
angle.  See  Figure  8.1,  where  point  B  is  replaced  by 
point  D. 

2.  No  return  can  be  made  to  points  which  have  just  been 
left;  see  Figure  8.2,  where  point  A  is  replaced  by 
point  D.  If  point  D  is  still  the  worst  point  of 
triangle  BCD,  then  point  B  Is  replaced  by  point  E. 

3.  If  the  best  vaiijed  vertex  remains  unchanged  for  more 
than  M  Iterations,  then  the  simplex  size  is  reduced 
by,  for  example,  halving  the  distance  of  all  other 
vertices  from  that  vertex.  The  next  stage  can  then 
sca.'t.  For  example,  as  in  Figure  8.3,  point  C  is 
sequentially  repeated  in  five  triangles  ABC,  BDC, 

DEC,  EFC,  and  FGC;  and  after  the  size  of  the  fifth 
triangle  FGC  is  reduced  to  triangle  F'G'C,  the  ru’'  I 
is  applied  at  the  new  triangel  F'G'C. 

The  magnitude  of  M  depends  on  the  number  of 
variables.  Spendley  [68]  suggested  that 


M  =  1.65n  +  0.05  n^ 


8.1 
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where  n  is  the  total  number  of  design  variables.  The 
search  can  finally  be  stopped  when  the  simplex  is  small 
enough  to  locate  the  optimum  adequately. 

It  should  be  noted  that  Rule  1  and  Rule  2  force  the 
simplexes  to  circle  continuously  about  an  indicated  opti¬ 
mum,  rather  than  oscillate  over  a  limited  range  such  as  a 
ridge  (see  Figure  8.2).  If  one  ensures  that  at  any  point 
violating  a  constraint,  a  sufficiently  large  negative 
value  (or  positive  value,  if  one  is  m^ocimizing)  is  set  for 
the  response  at  such  a  point,  the  system  of  simplexes  will 
then  move  along  rather  than  cross  the  constraints. 

The  advantages  of  this  method  are  as  follows; 

1.  It  is  easy  to  apply. 

2.  It  is  useful  when  analytical  or  numerical 
derivatives  of  the  objective  function  are  not 
available. 

The  disadvantage  of  this  method  is  that  it  could  not 
guarantee  that  the  global  optimum  design  is  obtained. 
Therefore,  several  different  initial  locations  (designs) 
should  be  considered  to  get  global  optimization. 


Energy  approach  and  laminated  plate  theory  approach 
are  both  utilized  in  damping  analysis.  The  main  advantag 
of  the  energy  approach  over  the  laminated  plate  theory 
approach  is  that  the  interlaminar  stresses  are  included, 
The  disadvantage  of  the  energy  approach  is  that  it  takes 
much  more  computer  time  than  the  laminated  plate  theory 
approach  does.  The  influence  of  interlaminar  stresses 
can  be  neglected,  if  the  thickness  of  each  group  of  lami¬ 
nated  composite  is  thin  (a  detailed  discussion  is  given 
in  Section  9.5).  Then  laminated  plate  theory  approach 
for  damping  analysis  is  a  good  way  for  the  optimization 
analysis  for  thin  structures. 
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Generally  speaking,  the  design  variables  of  continu¬ 
ous  or  discontinuous  fiber  composites  include  fiber  volume 
fraction  V£,  fiber  moduli  E^,  fiber  aspect  ratio  s/d,  the 
longitudinal  distance  between  fiber  tips  p,  and  fiber 
orientation  angle  6  for  certain  special  stacking  sequence. 
Among  those  six  variables,  only  fiber  aspect  ratio  s/d  and 
fiber  tip  longitudinal  distance  p  are  treated  as  the 
design  variables  in  the  optimization  analysis.  Since 
orthotropic  materials  are  utilized  In  the  optimization 
analysis,  the  fiber  direction  6  is  kept  as  certain  con¬ 
stants.  The  fiber  volume  fraction  is  also  set  as  a 
constant  (0.65)  jn  this  study,  since  the  fiber  composites 
are  made  from  prepreg  tape.  Stacking  sequence  is  fixed 
and  most  of  this  study  is  concentrated  on  graphite-epoxy, 
Kevlar-epoxy  and  the  hybrid  of  those  two  composite 
materials,  so  that  only  two  sets  of  fiber  moduli  are 
considered. 

Two  different  cases  are  considered  in  optimization 
analysis: 

Case  one  is  to  optimize  the  specific  stiffness  and 
material  damping  of  an  orthotropic  square  plate  simply 
supported  on  ail  four  sides  under  free  vibration. 

The  equation  of  motion  for  this  case  is 


8.2 
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where  the  Dj^j  are  derived  from  laminated  plate  theory,  and 
W  is  the  flexural  deformation.  Let 


W  =  -Nnn  «n  ^(t)  8.3 

By  separation  of  variables,  the  form  (8.4)  is  obtained  for 
T{t);  while  Fj^  (x)  and  (y)  satisfy  Equation  8.5 

T(t)  =  e  (i  t  k  =  1,  2,  3  .  .  .  . 

^  8.4 


''m  2(Di2  +  022 


8.5 


where  is  the  constant  to  be  determined,  and  i  is  J -1 
in  Equation  8.4. 

In  accordance  with  the  simple  supoort  bcundary 
conditions,  Fjjj  and  Kj^  functions  can  be  assumed  as 


sin 


miry 

a 


8.6 


sin 


nirx 

~15“ 


8.7 


where  a  and  b  are  the  lengths  of  the  edges  along  x  and  y 
directions,  respectively. 

By  substituting  Equations  8.6  and  8.7  into  Equation 
8.5,  one  obtains 


4  4  2,4  „4,4 

D^l  +  2(Di2  +  SDgg)  54^  +  D22 

a  abb 


8.8 


Letting  b=a  for  a  square  plate  and  considering  only 
the  first  mode  of  vibration  (i.e.,  ni-n=l),  one  can 
obtain  as 


T  ^*■-'11  ^^'12 


f  HU 
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The  natural  circular  frequency  of  the  first  mode,  w^,  is 
defined  as 


<0i 


8.10 


Then  the  first  mode  solution  of  Equation  8.2  is 
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If  the  viscoelastic  behavior  of  the  materials  is 
considered,  the  complex  form  of  should  be  included, 

i.  e. 


0)  = 


5  ‘“ti 


2D* 


4-  4D 
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D* 

^22' 


8.12 


After  separating  the  real  and  imaginary  parts  and 
neglecting  the  higher-order  terms  of  the  binomial 
expansion  of  the  quantity  on  the  right-hand  side  in 
Equation  8.12,  the  following  expression  is  obtained. 


*  .  .  « 

oj,  =(*),  +  i  a>, 

J.  X  X 


where 


(D|l  +  2D^2  +  "^^66  ^22^ 


8.13 


8.14 


^ 

Wl 

^  2a‘- 


d;i  .  2d;3  ^  4D;g  d’^ 

0“>ll  *  ♦  4D’5 


8.15 


I 


Hence,  for  viscoelastic  material.  Equation  8.11  should  be 


rewritten  as 
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.  /  4  TTX  .  ity  iuji  t,  a 

W  *  sin  ~  sin  ^  e  i  )  e  ■*•  8-io 

Equation  8.16  is  similar  to  the  mathematical  model  of 
logarithmic  decrement  {42],  and  the  logarithmic  decrement 
5  can  be  approximated  by 

8.17 

for  light  damping,  i.e.  when  u£/u)|  <  <1. 

By  using  Equation  2.22,  the  loss  factor  of  this 
system  is  then  determined  as 


2Jt 


Ui, 


w 


i- 


m: 


wi 


a  2ff 


(1) 


0) 


^il  '*~^Di2  '*^66  '*’  *^22 


8.18 


The  value  of  depends  on  the  stiffness.  High 

stiffness  reduces  the  deflection  at  resonance;  the 
material  damping  can  further  reduce  the  displacement  at 
resonance.  Gibson  (61]  shows  that  flexural  vibration  for 
a  double  cantilever  beam  under  forced  vibration  at 
resonant  frequency  depends  inversely  on  the  product  of 
material  damping,  area  moment  of  inertia  of  cross  section 
and  Young's  modulus.  It  should  be  noted  that  his  result 
is  based  on  the  equation  of  motion  for  free  vibration  and 
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the  boundary  conditions  of  forced  vibration.  For  the  case 
considered  here,  the  product  of  the  area  moment  of  inertia 
and  Young's  modulus  corresponds  to  the  generalized 
stiffness  D,  which  is  defined  as 

D*  =  D'  +  i  6"  8.19 

D'=  D'  +  20^2  ^  ^22 

B-  .  .  2D£2  *  4DJ5  .  8.21 

Equation  8.18  can  be  rewritten  as 

n  =  D"/D'  8.22 

In  order  to  have  small  resonant  deformation,  a  high 
value  of  the  product  of  stiffness  and  material  damping  is 
required.  But  high  specific  stiffness  is  the  major  pro¬ 
perty  of  fiber  composites  to  be  widely  used  in  space 
vehicle  and  aircraft.  Thoraforo,  the  mathematical  formu¬ 
lation  of  optimization  on  damping  and  specific  stiffness 
is  to  seek  the  maximum  value  of  the  objective  function 


fj^  (p/d,  s/d) 


8.23 


of  two  variables  p/d  and  s/d,  where  is  assumed  in  the  form 

0'  O' 


1 


(p/d,  s/d)  =  Tj  +  Tj 


*  n 


PO 


O‘o 

"po  *  Tti 


or 


8.24 


O'  PO 


f,  (p/^.  s/d)  = —  (T;  ♦  Tj -A) 


8.25 


The  ranges  for  design  variables  p/d  and  s/d  are 


0  <  o/d  <  0.05  s/d 


ft  ftft 


25  <  S/d  <  10000 


8.27 


Where  Tj  and  are  wieghting  constants,  is  the  damping  value  of 
epoxy,  p  is  the  density  of  designed  composite,  D  and  g 
and  n  of  design  composite  are  defined  In  equations  8,20  and  8.22, 
respectively,  and  Oq  and  po  are  the  corresponding  values  of  Djj  and 
density  of  unidirectional  continuous  graphite  reinforced  epoxy  having 
the  same  thickness  and  fiber  volume  fraction.  If  high  specific  stiffnes 
is  important  in  structures,  could  be  chosen  as  a  high  value. 
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Alternatively,  if  small  resonant  deformation  is  the  major 
consideration,  T2  should  be  a  high  value. 

Case  two  is  to  optimize  the  specific  stiffness  and 
material  damping  of  an  orthotropic  square  plate  clamped  on 
all  four  sides  under  free  vibration. 

The  equation  of  motion  is  just  the  same  as  Equation 
8.2.  The  first  mode  solution  is  assumed,  according 
to  reference  (701 

K  =  F(x)  H(y)  T(t)  8.28 

where 


F(x)  =  cos  X^x  -  cosh  Xj^x  +  sin  Xj^x  -  sinh 

8.29 


H(y)  =  Bj^  cos  Xj^y  -  Bj^  cosh  Xj^y  +  sin  Xj^y  -  sinh  Xj^y 

8.30 


sin  Xj^a  -  sinh  Xj^a 
-  cos  Xj^a  +  cosh  Xj^a 


8.31 


and  Xj^  is  the  constant  to  be  determined. 

The  natural  frequency  of  first  mode  vibration  of  a 
clamped  square  plate  was  given  in  reference  (70)  as 


i  (5.14Dj^^  +  1.55(2022  +  40^^)  +  5. 140221** 


8.32 
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Ey  a  similar  approach  to  that  prescribed  in  case  one,  one 
will  obtain  the  loss  factor  of  the  system  for  case  two  as 


+  1.55  (2DJ2+  4D^g)  +  5.140^2 

3.14b|"”+‘  1.55  '(2d”|2+  +  5. 140^2 


The  mathematical  formulation  of  optimization  for  this 
case  is  to  seek  the  maximum  value  of  the  objective 
function 


£2  (p/d,  s/d)  8.34 

of  two  variables  p/d  and  s/d,  where  £2  defined  in 
Equation  8.37 

The  ranges  of  design  variables  are 


0  S  p/d  S  0.05  s/d  8.35 

25  S  s/d  5  10000  8.36 

where 


f2  (Pi*  s/d) 


8.37 


and 


CHAPTER  9 

RESULTS  AND  CONCLUSIONS 


9.1  Preliminary  Remarks 

The  damping  analyses  presented  in  Chapters  3  to  6  are 
applicable  to  all  kinds  of  fiber  composites,  including 
continuous  fiber  composites,  discontinuous  fiber  compo¬ 
sites,  symmetrically  or  unsymmetrical ly  laminated  or  lami¬ 
nar  composites,  randomly  oriented  fiber  composites,  etc. 
The  optimization  analysis  on  damping  and  specific  stiff¬ 
ness  presented  in  Chapter  8  is  based  on  orthotropic 
material.  However,  a  similar  approach  is  applicable  for 
more  general  anisotropic  materials.  Since  most  widely- 
and  practically-used  fiber  composites  are  symmetrically 
laminated  fiber  composites,  the  numerical  analysis  of 
this  study  is  concentrated  on  certain  kinds  of  symmetri¬ 
cally  laminated  fiber  composites,  unidirectional  laminar 
composites,  and  in-plane  randomly  oriented  short-fiber 
composites. 

As  mentioned  before  (see  Chapters  3  and  5),  damping 
njj  is  defined  as  (or  .  Since  both  e;^  and  e;; 

(or  and  D^)  are  functions  of  E^,  s/d,  0,  P  (the 

longitudinal  distance  between  fiber  tips),  Vj  and  etc., 
the  variations  of  rij£  and  E^  (or  D^)  may  not  follow  the 
same  pattern.  Consequently,  in  the  numerical  results. 
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three  dependent  variables  (i.e.,  n'y.,  E'^^,  and  E^)  are  dis¬ 
cussed  in  laminar  composites,  while  more  dependent 
variables  are  involved  in  other  composites.  Among  the 
independent  variables,  most  interest  is  focused  on  E£,  Ej^, 
s/d,  and  0.  P  is  treated  as  an  independe.it  variable  only 
for  optimization,  and  is  set  as  a  constant  in  most 
cases . 

The  numerical  results  presented  in  this  chapter  are 
based  on  the  following  assumptions: 

1.  The  matrix  epoxy  behaves  as  a  linear,  isotropic, 
viscoelastic  material  [15]. 

2.  All  fibers  are  also  linear  viscoelastic  materials,  and 
damping  of  fibers  is  independent  of  direction. 

3.  The  dependence  of  modulus  on  frequency  is  not 
considered  at  the  present  time. 

9.2  Damping  and  Stiffness  of  Unidirectional 
Fiber  Composites 

The  numerical  results  of  damping  and  stiffness  analy¬ 
sis  presented  in  Chapter  3  are  discussed  in  this  section. 
It  should  be  noted  that  the  modified  Cox's  model 
(i.e.  Equations  3.8  and  3.13)  is  used  to  calculate  the 
longitudinal  modulus  of  composite  materials. 

Figures  9.1,  9.2,  and  9.3  present  the  nondimensional 
plots  of  and  nx/h^  cf  graphite  epoxy 

composite  as  functions  of  fiber  aspect  ratio  s/d,  using  0 
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(the  angle  between  fiber  direction  and  loading  direction) 
as  paraiiieter.  For  small  angle,  say  Q<10*,  the  storage 
modulus  sharply  increases  as  fiber  aspect  ratio  increases; 
the  loss  modulus  reaches  its  maximum  around  s/d=130  and 
then  reduces  as  s/d  increases,  while  damping  apparently 
decreases  as  s/d  increases.  Comparing  0®  curves  with  the 
30®  curves  in  those  drawings,  one  can  notice  that  when 
0=0®  low  damping  at  large  s/d  is  caused  mainly  by  high 
storage  modulus.  Comparing  the  60*  curves  with  the  90' 
curves  in  those  drawings,  it  is  recognized  that  for  large 
angles,  the  differences  in  loss  moduli  have  more  influence 
on  damping.  Alternatively,  Figures  9.4,  9.5,  and  9.6  show 
the  nondimensional  plots  of  /H^’^  ,  /E'^  ,  and  hj^/hn,  of 

graphite  epoxy  composites  as  functions  of  0,  using  s/d  as 
parameter.  It  is  clearly  observed  that  tl»^  shorter  the 
fiber,  the  stronger  is  the  dependence  of  E^/EjJj  ,  ' 

and  rijj/Ojjj  on  s/d  for  small  angles,  say  0<1O.  When  0^45, 
there  is  not  any  change  of  ,  EJ^/e'^  ,  and  nx/n^i 

caused  by  changing  s/d.  In  addition,  Figure  9.5  shows 
that  the  maximum  value  of  occurred  between  0=C® 

and  0=15®  for  fiber  aspect  ratio  ranging  from  25  to  10000„ 
Similar  properties  of  unidirectional  boron  epoxy 
composite  and  glass  epoxy  composites  are  observed.  But, 
for  0=0,  the  maximum  values  of  E^/El^  occur  around  s/d=500 
and  s/d=50  for  boron  epoxy  composite  and  glass  epoxy 
composite,  resp ->ctively.  Also,  the  damping  of  boron  epoxy 
composite  is,  in  general,  higher  than  that  of  graphite 
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epoxy  composite.  This  is  because  boron  itself  has  higher 
dcunping  than  graphite  does. 

The  corresponding  drawings  of  Figures  9.1,  9.2,  and  9.3 
for  unidirectional  Kevlar  epoxy  composite  are  given  in 
Figures  9.7,  9.8,  and  9.9,  respectively.  It  ii  observed 

that  Kevlar  epoxy  composite  has  the  highest  material  damp¬ 
ing  among  those  four  kinds  of  composite.  This  is  due  to 
the  fact  that  Kevlar  (aramid)  is  a  kind  of  polymer  which 
usv.>lly  has  high  damping.  It  is  interesting  to  notice 
that  of  Kevlar  epoxy  composite  monotonical ly 

increases  with  respect  to  s/d  when  0=0  “.  All  other  proper¬ 
ties  of  Kevlar  epoxy  composite  are  similar  to  those  of 
graphite  epoxy  composite. 

Figure  9.10  presents  the  plots  of 
of  Kevlar  epoxy  composite  as  function.s  of  0  with 
s/d=l00.  This  figure  shows  that  high  damping  of  such 
material  is  induced  by  the  fact  that  the  values  of 
and  /E^  are  very  close.  The  damping  of  Kevlar  epoxy 
composite  Is  insensitive  to  angle  change  because  the 
reduction  rates  of  and  are  very  close.  Figure 

9.11  shows  the  plots  of  and  E^/E’^  ,  and  of 

Kevlar  epoxy  composite  and  graphite  eooxv  composite  as 
functions  of  0,  while  s/d  is  kept  as  lOU.  This  suggests 
idea  that  if  a  hybrid  composite  (graphite  epoxy  and  Kevlar 
epoxy  composite)  is  made,  damping  of  this  hybrid  composite 
will  be  higher  than  that  of  graphite  composite,  while  the 
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stiffness  of  the  hybrid  composite  will  be  greater  than  that 
of  Kevlar  epoxy  composite. 

Figure  9.12  presents  of  graphite  epoxy  composite 

as  function  of  0  for  different  fiber  volume  fractions 
when  s/d=100.  Differences  in  caused  by  different 

volume  fractions  decreases  when  angle  increases.  It  is 
also  observed  that  an  apparent  reduction  in  occurs 

when  angle  increases  from  0®  to  15®.  Figure  9.13  exhibits 
the  plots  of  for  graphite  opoxy  composite  versus  0 

for  different  fiber  volume  fractions,  when  s/d=100.  The 
peak  value  of  damping  moves  from  about  50®  to  30®  when  the 
volume  fraction  changes  from  0.1  to  0.6. 

Three-dimensional  plots  using  and  as  vertical 
axis,  and  s/d  and  0  as  two  base  axes  of  Kevlar  epoxy 
composites  are  given  in  Figures  9.14  and  9.15, 
respectively.  The  associated  contour  curves  are  presented 
in  Figures  9.16  and  9.17. 

9.3  Damping  and  Stiffness  of  Randomly  Oriented 
Short-Fiber  Composites 

Figure  9.18  gives  the  nondimensional  plots  of 
E^E'j^  ,  and  E"j./E"jjjOf  glass  epoxy  composite  as  functions  of 
s/d.  It  is  observed  that  damping  reaches  its  maximum  value 
at  small  values  of  s/d  and  decreases  as  s/d  increases.  The 
storage  modulus  E^/E^  monotonical  ly  Increases  as  s/d 
increases.  However,  loss  modulus  almost  remains 

unchanged.  In  Figures  9.19  and  9.20,  nondimensional  plots 
of  ^r/^n^,  E'j,/E^  and  E^EJ^  for  graphite  epoxy  composite 


and  Kevlar  epoxy  composite  are  presented  separately.  The 
variations  of  nj-Zn^.  and  for  Kevlar  epoxy  composite 

or  graphite  epoxy  composite  are  similar  to  those  of  glass 
epoxy  composite.  But,  for  Kevlar  epoxy  composite,  E^/F^ 
increases  as  s/d  increases,  while  graphite  epoxy 
composite,  E"/E^  has  its  maximum  value  around  s/d=120. 
Figures  9.21,  9.22  and  9.23  show  similar  plots  of  Hcr/^Gm 
,  and  G^/G^  as  functions  of  s/d  for  glass  epoxy 
composite,  graphite  epoxy  composite,  and  Kevlar  epoxy 
composite,  respectively.  Variations  of  Pcr/’^Gm'  ^’r^^ 
and  G^/G^  are  almost  the  same  as  those  of  nj-ZOm,  ^r/^ra' 
and  E“/E^  of  corresponding  composite,  respectively. 

Figure  9.24  shows  the  nondimensional  plots  of  ny/ 
as  function  of  fiber  volume  fraction  Vf,  using  s/d  as 
parameter.  It  is  observed  that  nj./njj,  decreases  as 
increases  for  small  s/d  as  well  as  for  large  s/d. 

Three-dimensional  plots,  using  Ev  and  as  vertical 
axis,  and  Ef  and  s/d  as  two  base  axes,  of  four  kinds  of 
fiber  composite  are  given  in  Figures  9.25  in  9.26,  sepa¬ 
rately,  The  associated  contour  curves  are  presented  in 
Figures  9.27  and  9.28. 

9.4  Damping  and  Stiffness  of  Laminated  Composite — 
Laminated"  P~1  ate  Theory  Approach 

The  numerical  results  of  the  damping  and  stiffness 
analysis  described  in  Chapter  5  are  presented  in  this 
section.  Five  types  of  lamination  (quasi-isotropic 
lamination,  angle  ply  lamination,  cross  ply  lamination. 
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unidirectional  lamination,  and  a  special  kind  of  lamination 
g)  of  graphite  Epoxy  composites  are  considered  in 
numerical  examples. 

Figure  9.29  exhibits  nondimensional  plots  of  Dii  and 
F^ll  normal  damping)  for  [02/902/452/'*452 1  g, 

(452/02/-452/902ls»  {452/'*4  2/502^®2^s  functions  of 

s/d.  It  should  be  noted  that  all  those  three  kinds  of 
lamination  have  the  same  values  of  in-plane  longitudinal 
stiffness  and  in-plane  longitudinal  damping  jnn*  But 
Figure  9.29  shows,  as  expected,  that  for  the  same  s/d, 
[O2/9O2/452/-452 1  g  has  the  highest  value  of  and  the 
lowest  value  of  f452/-452/902/02  ]  s  lowest 

value  of  and  the  highest  value  of  pn^,  while  [452/O2/ 
-452/90213  has  almost  the  average  values  of  and  of 

of  the  other  two  kinds  of  lamination.  Alternatively, 
Figure  9.30  shows  the  nondimensional  plots  of  and  phgg 
(twisting  damping)  for  those  three  kinds  of  lamination  as 
functions  of  s/d.  In  this  drawing  for  same  s/d,  [O2/9O2/ 
452/-452]g  has  the  lowest  value  of  E^g  and  the  highest 
value  of  pHgg/  5  highest  value  of 

Dgg  and  the  lowest  value  of  ph5g,  and  [452/02/“452/902] g 
still  remains  around  the  average  values  of  those  other  two 
kinds  of  lamination.  Figures  9.29  and  9.30  actually  show 
us  how  to  design  the  quasi-isotropic  laminated  plate  for 
either  maximum  damping  or  maximum  stiffness  purpose. 

Figures  9.31  and  9.32  show  the  nondimensional  plots  of 
and  (in-plane  longitudinal  dcimping)  of  angle  ply 
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lamin-^ted  graphite  epoxy  composite  as  function  of  G,  using 
s/d  as  a  parameter.  increases  with  s/d  for  small 

angles,  say  0<2O*.  The  variation  of  with  s/d  reduces 
as  0  increases,  >.  nd  very  small  differences  exist  for  G>60®. 
The  distinctively  decreases  as  s/d  increases  for 

small  angles  (say.  e<30“),  while  the  variation  of 
with  s/d  reduces  as  0  Increases.  Figures  9.31  and  9.32 
are  comparable  to  Figures  9.4  and  9.6  correspondingly-  The 
reduction  in  jnn  (or  the  increment  in  of  angle  ply 

laminated  composite  is  not  as  sharp  as  that  of  unidirec¬ 
tional  composite,  when  s/d  is  increased  and  the  angle  is 
small . 


In-plane  shear  stiffness  and  in-plane  shear 
damping  j-rigg  of  angle  ply  graphite  epoxy  composite  as 
function  of  6,  using  s/d  as  parameter  are  shown  in  Figures 
9.33  and  9.34,  separately.  Maximum  and  minimum  jTigg 

occur  at  0=45'.  Larger  dependences  of  ^66  and  jHgg  on  0 
are  observed  for  larger  s/d.  The  reason  why  /D'^/ 

Dee/^Gm'  F’lll/'^m'  presented  in 

drawings  is  that  those  values  are  very  close  to  the  values 
of  E^^/e;  ,  E^./g;  ,  inu/Hni.  and  jneS^^Cm  separately. 

Figures  9.35-9.38  present  the  nondimensional  plots  of 


F^ll'  ^66'  F^66  graphite  epoxy  composite 

laminated  with  four  kinds  of  lamination  ([45/-45]4g, 
I452/O2/-452/9O2IS/  (0/90]4g,  and  lOglg)  as  functions  of 

s/d.  It  is  observed  that  for  the  same  s/d,  [45/~45]4g  has 
the  highest  values  of  D^g  and  and  the  lowest  values 


mmmm 


V  V  V  V  v 


93 


of  and  ftigg.  The  values  of  and  pn65  for  (0/90l4g, 
and  for  [Oglg  same,  and  those  values  are  almost 
Independent  of  s/d.  This  is  because  E,p  and  of 

fiber  composites  are  assumed  to  be  indendpent  of  the  fiber 
length.  D  1^2^  of  (Og]g  is,  as  expected,  higher  than  of 
(0/90l^g  at  the  seune  s/d,  but  phii  of  those  two  kinds  of 
lamination  are  surprisingly  close.  The  ,  pn_L]^, 
and  phgg  of  [452/02/-452/902l4g  are  just  somewhere  near 
the  middle  range  of  those  of  [45/-45]4g  [0/90] 

The  characteristic  property  of  [0n'^“®4n^s  that 
its  and  are  very  close  to  zero.  For  example, 

D^g/D'ij^  or  0*25/0*2^2^  of  this  lamination  for  graphite  epoxy 
composite  is  less  than  0.02  for  all  angles,  while  those  of 
angle  ply  lamination  are  less  than  about  0.11  for  all 
angles.  Numerical  results  show  that  D{2.»  f'^11'  *^66' 
pfigg  of  [®2/~®8^s  very  close  (difference  within  2%)  to 
those  of  [GZ-Qlgg.  But  Ej^2.  ^nd  E55  (or  jn2.1 
are  in  general  less  (or  greater)  than  those  of  [©/-Olgg 
(maximum  differences  ranging  from  13%  to  30%) . 

9.5  Damping  and  Stiffness  of  Laminated  Fiber  Composite — 

Energy  Approach 

There  are  twenty-one  different  kinds  of  material 
damping  induced  from  C^j /C^ j  (i,  j  =  1,  2  ...  6)  for 
anisotropic  materials.  In  the  energy  approach,  material 
damping  of  structure  is  just  the  ratio  of  the  sum  of 
^  c^j  ^ij  ®f  ^  ^ij  ^ij-  Therefore,  different 

values  of  material  damping  could  be  observed  for  one 


material  under  different  kinds  of  loading  and  boundary 
conditions.  In  the  author’s  opinion,  damping  analysis 
through  the  energy  approach  is  more  acceptable  than  lami¬ 
nated  plate  approach,  because  the  energy  approach  consi¬ 
ders  the  influence  of  all  stresses  on  the  dissipated 
energy,  while  laminated  plate  approach  considers  one 
stress  at  a  time.  The  results  of  damping  analysis  through 
the  energy  approach  presented  in  Chapter  6  are  discussed 
in  this  section. 

Figure  9,39  shows  the  difference  between  material 
damping  through  the  energy  approach  and  through  the  lami¬ 
nated  plate  approach  for  a  60x60x10  cross  plied 
(05/905/05)  graphite  epoxy  composite  under  uniform  in¬ 
plane  tension  load  along  x  direction.  The  difference  is 
mainly  caused  by  the  fact  that  there  exists  some  stress 
other  than  for  example,  negligible  in 

certain  regions.  The  difference  increases  when  the  value 
of  s/p  is  increased.  In  other  words,  continuous  fiber 
composites  have  more  apparent  difference  of  damping  values 
through  those  two  approaches  than  discontinuous  (short- 
fiber)  composites  do. 

Figure  9.40  presents  the  influence  of  plate  thick¬ 
ness  and  of  the  lamination  on  the  damping  value  by  the 
energy  approach  under  in-plane  uniform  tension  load.  It  is 
noticed  that  either  increasing  the  number  of  interfaces  or 
reducing  the  total  plate  thickness  will,  in  general,  reduce 
the  difference  between  damping  values  of  those  two 
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approaches.  It  should  also  be  noted  that,  for  an  isotro¬ 
pic  material  under  tensional  load  along  x  direction,  the 
damping  value  of  the  energy  approach  is  same  as  the  value 
of  through  the  laminated  plate  theory  app;:oach. 

More  numerical  results  of  using  the  energy  approach 
are  under  consideration.  In  particular,  the  effects  of 
Interlaminar  stresses  on  damping  and  stiffness  of 
laminated  composites  due  to  the  difference  in  stacking 
sequence,  mismatch  at  the  interfaces,  and  width-to- 
thickness  ratio  a/h  will  be  studied. 

In  addition,  we  are  also  interested  in  studying  the 
effects  on  damping  due  to  various  boundary  conditions 
under  the  same  flexural  load.  The  results  will  be 
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9.6  Experimental  Results  o f  Damping 
and  Stiffness 

Figure  9.41  presents  some  experimental  results  by  the 
impulse  hammer  technique  [35,  36]  about  material  damping 
and  stiffness  of  graphite  epoxy  composite  beams.  Experi¬ 
mental  results  of  stiffness  are  very  close  to  those  of  the 
analytical  results,  while  experimental  results  of  damping 
are  not  so  close  to  the  analytical  results.  A  similar 
trend  is  also  presented  in  Figure  9.42  where  experimental 
results  [35,  36]  about  damping  and  stiffness  of  off-axis 
unidi rectional  continuous  graphite  epoxy  composite  beams 
are  compared  with  analytical  results.  Tabic  9.1  shows  the 
experimental  results  of  damping  and  of  the  first  mode 
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nature  frequency  (thty  quantity  is  proportional  to  the 
stiffness  of  the  structure)  of  five  cross  ply  short-fiber 
glass  epoxy  cciiiposite  plates.  It  is  observed  that  high 
value  of  s/p  (either  increase  fiber  length,  s,  or  reduce 
the  gap  between  fiber  tips,  p)  will  increase  the  nature 
frequency  and  decrease  the  damping. 

9.7  Optimization  of  Deunping  and  Specific  Stiffness 

of  Fiber  Composites 

According  to  the  prescription  about  optimizing  damp¬ 
ing  in  Chapter  8,  some  optimum  designs  are  discussed  in 
this  section  for  graphite  epoxy  composites  and  for  hybrid 
composites  keeping  equal  to  0.65.  The  design  vari¬ 
ables  include  the  fiber  aspect  ratio,  s/d,  and  the  ratio 
of  the  longitudinal  distance  between  fiber  tips,  p,  to  the 
fiber  diameter,  d. 

Table  9.2  shows  the  influence  of  weighting  constants, 
Tj^  a. id  T2,  on  the  optimum  design  for  0/90/0  graphite  epoxy 
composites  under  situations  of  case  one.  By  increasing 
T2,  s/p  is  reduced;  while  by  increasing  Tj^,  s/p  is 
increased.  These  imply  that  to  have  smaller  resonant 
deformation,  larger  distance  between  fiber  tips  is 
preferred;  while  to  have  higher  specific  stiffness, 
smaller  distance  between  fiber  tips  is  required. 

Table  9.3  shows  the  optimum  design  for  cross  ply 
graphite  epoxy  composites  and  hybrid  composites  under  case 
one  situation  and  having  the  same  total  thickness,  with 
T]^=l  and  T2-5.  There  is  no  difference  in  the  optimum 
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values  of  design  variables  s/d  and  p/d  for  0/90/0  and 
0/90/0/90/0  graphite  epoxy  composites.  For  the  same  total 
thickness,  with  Tj^=l  and  T2=5,  discontinuous  fibers  give  a 
high  value  of  objective  function  in  the  stacking  sequence 
0*^/90^/0^/90^/s;  while  continuous  fibers  give  high  value 
of  the  objective  function  in  the  stacking  sequence 
oK/9oK/oG/9oG/9oG/g  where  the  superscripts  G  and  K  repre¬ 
sent  graphite  and  Kevlar,  respectively,  and  the  fiber 
composites  having  s/p=2593.  can  be  treated  as  continuous 
fiber  reinforced  composites.  When  comparing  the  value  of 
objective  function  of  each  optimum  design  for  each  case, 
one  can  discover  that  0^/90^/0^/90^/s  has  the  highest 
value  of  objective  function.  This  indicates  that  the  con¬ 
tinuous  hybrid  fiber  composite  plate  with  the  stacking 
sequence  0*^/90^/0^/90^/s  is  the  optimum  design  for  case 
one. 

Figure  9.43  shows  the  contour  curves  of  objective 
function  of  case  one  vs  s/d  and  s/p  (keeping  Tj^=l,  T2=5). 
It  is  observed  that  there  probably  exists  an  acceptable 
region  as  optimum  design.  In  fact,  the  term  (s/p)^^^  in 
TabJes  9.2,  9.3  and  9.4  is  the  average  value  of  ratio  s  to 
p  at  all  points  whose  objective  value  is  within  99.95%  and 
100.05%  times  of  the  objective  value  of  the  optimum  design 
shown  in  s/d  and  p/d  items.  Numerical  results  show  this 
region  ranging  from  s/d=2000  to  s/d  over  9000  for  case 
one.  That  means,  as  long  as  s/p  is  kept  equal  to  53.1  and 
s/d  over  2000,  one  will  get  an  acceptable  optimum  design 
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(when  Tj^al  and  T2=5)  for  0/90/0  or  0/90/0/90/0  cross  ply 
graphite  epoxy  composite  plates  under  case  one  conditions. 

The  optimum  designs  of  case  two  (see  Table  9.3)  is 
very  close  to  those  of  case  one,  although  D  for  those  two 
cases  are  different.  It  is  because  02^2 
small  relative  to  and  D22  (see  Equation  8.20  and 
Equation  8.38) . 


9.8  Conclusions 

Since  the  material  damping  and  stiffness  of  continuous 
or  discontinuous  fiber  composites  depends  upon  many  quanti¬ 
ties  such  as  stiffness  ratio  E^/Ejj^,  fiber  aspect  ratio 
s/d,  longitudinal  distance  beetween  fiber  tips  p,  loading 
(or  fiber)  angle  6,  fiber  volume  fraction  V^:,  fiber  and 
matrix  damping  values,  and  staclcing  sequence,  it  is  not 
feasible  to  discuss  the  material  damping  and  stiffness  by 
considering  all  the  possible  variations  of  all  the 
independent  quantities  simultaneously.  However,  based 
upon  the  numerical  results  presented  in  this  study,  we  can 
make  the  following  conclusions: 

A.  About  Unidirectional  Fiber  Composites 

1.  High  damping  could  be  achieved  by  small  0,  small 
s/d,  and  high  p;  while  high  stiffness  could  be 
obtained  through  high  Vj,  small  0,  large 

s/d,  and  low  p. 

2.  Damping  of  discontinuous  fiber  composites  is  no 
less  than  that  of  continuous  fiber  composites. 
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3.  For  small  angle,  high  s/d  gives  high  stiffness; 
while  low  s/d  brings  high  damping.  For  large 
angles,  damping  and  stiffness  are  independent  of 
s/d. 

B.  About  In-Plane  Randomly  Oriented  Fiber  Composites 

1.  Relative  to  unidirectional  fiber  composites,  in¬ 
plane  randomly  oriented  fiber  composites  have 
lower  stiffness  at  the  price  of  keeping  high 
damping  in  all  directions  parallel  to  the  plane. 

2.  Variations  of  shear  damping  and  shear  modulus  vs 
s/d  are  very  close  to  those  of  in-plane  damping 
and  in-plane  stiffness,  respectively. 

3.  Influence  of  on  In-plane  damping  becomes  less 
when  s/d  is  smaller. 

C.  About  Laminated  Fiber  Composites 

1.  In  a  laminated  composite,  there  are  eighteen 
damping  coefficients  in  general.  For  symmetric 
laminates,  the  number  of  damping  coefficients  is 
reduced  to  twelve.  Among  the  twelve,  there  are 
only  four  important  coefficients  namely; 

(or  t^22^  (extensional ) ,  ih65  (in-plane  sheer), 
pHxx  (flexural  normal),  and  (flexural  shear 

or  twisting)  . 

2.  As  unidirectional  composites,  laminated  composites 
with  higher  damping  are  usually  with  lower 


stiffness . 
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3.  In  general,  laminates  with  large  fiber  aspect  ratio 
s/d  and  very  small  fiber  tip  spacing  (p=0)  will 
produce  high  damping  and  high  stiffness. 

4.  The  analytical  results  in  damping  from  the  energy 
approach  are  higher  than  the  corresponding  results 
obtained  from  laminated  plate  approach.  This 
implies  that  the  interlaminar  stresses  tend  to 
dissipate  more  energy. 

5.  The  influence  of  interlaminar  stresses  on  damping 
is  very  significant  for  the  laminates  with 

large  fiber  aspect  ratio  s/d  or  relatively  large 
thickness- to-width  ratio  h/a. 

RECOMMENDATI ONS 

1.  Use  aligned  fiber  with  large  fiber  aspect  ratio 
s/d  and  very  small  fiber  tip  spacing  (p=0) 

(if  fabrication  process  permits) . 

2.  Use  hybrid  fiber  composites.  In  particular,  use 
the  conibination  of  Kevlar-epoxy  and  graphite- 
epoxy  laminated  composites.  The  advantage  of  using 
such  a  combination  is  due  to  the  higher  damping  in 
Kevlar  fibers  and  higher  stiffness  in  graphite 
fibers . 

About  Optimization 

Optimumization  of  damping  and  specific  stiffness 

depends  on  the  preference  of  high  specific  stiffness 

or  high  damping  (Tj^  or  T2),  material  property  of 

fiber,  the  stacking  sequence  for  hybrid  composite. 
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and  boundary  conditions.  A  continuous  fiber  rein¬ 
forced  hybrid  orthotropic  composite  plate  with  stack¬ 
ing  sequence  0^/90^/0^/90^  is  the  optimum  design  for 
clamped  (or  simply  supported)  square  plate  under  the 
free  vibration. 
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Figure  9.12:  Plots  of  E’/E'  vs  0  Using  V£  as  a  Parameter  for  Graphite  Epoxy  Coroposi 
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Table  9.1:  Experimental  Results  of  Og/90^/0^  Glass 
Epoxy  Composite  Plates 


s  (mm) 

p(mm) 

cl(mm) 

Vf 

f^.(Hz) 

n 

19. 

1.59 

.009 

.558 

92.8 

.0042 

25.4 

1.59 

.009 

.565 

95.3 

.0041 

38.1 

1.59 

.009 

.575 

97.1 

.0040 

25.4 

0.05 

.009 

.579 

109.8 

.0036 

Table  9.2:  Influence  of  Weighting  Constants  on  Optimum 
Design  for  Case  One 


Lamination 

•^1 

^2 

s/d 

pTd 

^^ave 

D’po 

J1 

pDi 

0/90/0 

1 

3 

9022. 

121. 

74.4 

1.27 

.47 

0/9010 

1 

5 

9022. 

170. 

53.1 

1.19 

.52 

0/90/0 

1 

10 

9008. 

214. 

42.1 

1.12 

.56 

0/90/0 

0 

5 

9020. 

267. 

33.8 

1.06 

.60 
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Table  9.3:  Optimum  Designs  of  Cross  Ply  Composite 
Plates  for  Case  One 


Lamination 

la 

■~s7a — 

av'e 

D'p. 

n 

^'m 

pDi 

(0/90/0)^ 

1 

5 

9022. 

170. 

HHI 

1.19 

.52 

(0/90/0/90/0)^ 

1 

5 

8994. 

169. 

53.1 

1.19 

.52 

0^/90^Vo’^/90^/s 

1 

5 

9007. 

152. 

59.4 

1.29 

.53 

0*^/90^/0‘^/90^/s 

1 

5 

9993. 

3.2 

2593.1 

.69 

Table  9.4:  Optimum  Designs  of  Cross  Ply  Composite 
Plates  for  Case  Two 


Lamination 

'i'l 

*2 

s/d 

p/d 

(s/p) gyg 

D'p* 

H 

pDS 

^m 

(0/90/0)^ 

g 

9010. 

170. 

53.0 

5.10 

.46 

(0/90/0/90/0)^ 

g 

9020. 

170. 

53.0 

5.10 

.•46 

0^/90^7oV90*^/s 

1 

5 

9010. 

152. 

59.3 

5.53 

.48 

0*'/90^/0®/90‘^Vs 

1 

5 

9993. 

3.2 

2747.8 

4.80 

.65 

K  /  K  ,_G  ,-_G  , 
/90  /O  /90  /s 


2747.8  4.80 
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APPENDIX  A 

ENERGY  EXPRESSION  OF  DAMPING  FOR  LAMINATED  COMPOSITES 


Equation  2.22  indicates  that  the  material  damping,  ri/ 
equals  to  the  ratio  of  the  energy  dissipated  in  a  cycle, 
(UjjJcyc^  the  product  of  2^,  and  the  maximum  strain  energy 
stored  in  the  body,  Ug,  under  a  periodic  vibration.  Since 
the  maximum  strain  energy  of  a  beam  under  static  load  and 
under  first  mode  vibration  are  equal,  if  the  maximum 
deflection  are  the  same.  If  the  material  damping  of  this 


beam  id  known. 


tite  energy  uxssxpateu  in  the  ueam  within 


cycle  of  the  first  mode  vibration  can  be  determined  by 


^^D^cyc  =  2ti  n  Ug  A.l 

where  Ug  can  be  determined  by  calculating  the  strain 
energy  of  the  statiscal ly-deformed  beam  having  the  same 
amplitude  of  the  deflection  as  that  of  the  first  mode 
vibration.  As  it  is  well  known  that  the  dominant  stress 
in  a  beam  under  flexure  vibration  is  normal  stress,  say 
Therefore,  Ug  and  can  be  obtained 

approximately  as 


14  6 


and 


14  7 


eye 


2a 


°x 


A. 3 


Where  Gj^  is  the  normal  strain,  and  Is  determined  by 
E^/E^  .  To  be  more  accurate,  all  stresses  and  strains  in 
the  beam  should  be  included  when  calculating  the  maximum 
strain  energy  and  the  dissipated  energy. 

Ungar  and  Kerwin  [56]  derived  the  expressions  for  the 
loss  factor  of  any  series-parallel  array  of  m  viscoelastic 
springs 


n 


m 

^^Dj^cyc 

3°1 

ra 

S  2tt 


A. 4 


A. 5 


where  U  ^  denotes  the  maximum  strain  energy  stored  in  a 
deflection  Xj  in  the  jth  lossless  spring  of  stiffness  kj, 
and  denotes  the  energy  dissipated  in  the  jth 

viscoelastic  spring  during  a  cycle.  Since  there  are  six 
different  components  in  the  stiffness  matrix  for  an 
orthotropic  material,  one  can  treat  an  orthotropic 
material  as  a  structure  composed  of  six  springs. 
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Consequently,  the  material  damping  of  a  laminated  n  plies 
composite,  n»  in  the  first  mode  vibration  can  be 
determined  by 


In  which  Vg  is  the  volume  of  the  qth  layer  in  the 
laminated  composite,  <^j  ,  C^j  and  are  the  storage 
moduli,  the  loss  moduli,  and  the  complex  moduli  for  a 
viscoelastic  material,  respectively- 


APPENDIX  B 

FORMULATION  OF  FINITE  ELEMENT  METHOD 

This  appendix  briefly  prescribes  the  formulation  of 
the  finite  element  method  utilized  in  energy  approach  to 
determine  the  displacement  and  the  strain  fields  for  an 
elastic  body.  The  material  constants  are  also  presented 
here. 


B.l  Stiffness  Matrix 

The  nonzero  elements  of  stiffness  matrix  [C]  for  an 
orthotropic  material  when  its  coordinates  coincide  with 
principal  coordinates  are  given  in  Reference  (57]  as: 

—  (1  ^23  ^ 

^12  "  ^''12  ^^32  ’^13^®2  ^ 

^13  "  ^"^13  '*■  ''.12  ''23^%  ^ 

C22  =  (1  -  Vj^3  V33^)E2  a  B.4 

^23  “  ^''23  '’21  ''13^^3  ^ 

C33  =  (1  -  V3^2  ''21^^3  ^ 
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15  0 


^44  ~  ^23 


B.7 


C55  ~  *^31 


B.8 


^66  “  *^12 


B.9 


where 


~  ^12^21  “  ^23^32  “  ^31^X3  *■  ^'^2l''32''l3  b.10 

Since  the  material  properties  of  a  laminar  fiber  composite 
with  the  coordinates  parallel  to  and  transverse  to  the 
fiber,  as  Figure  B.l,  can  be  treated  as  transversely 
(in  2,  3  directions)  isotropic. 


Figure  B.l;  Principal  coordinates  of  fiber  composite 
material 
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one  has 

El  =  El 

^2  “  E>p 

E3  = 

B.ll 

^12  ^  <^LT 
^11  "  ^LT 


'23 


"2  ( 1  +  v,p^  I ) 


and 


^12  "  ^LT 


V21 


V 


LT 


E 


T 


V23  -  V.j.^1 


^^32  “  '^TT' 


^13  -  ^LT 


''31“  '’TL 


B.12 


I 


15  2 

where  ^LT'  ^IT  defined  by  Equations  3.13, 

3.15-3.17,  and  is  assumed  to  be  the  Poisson's  ratio 

of  the  matrix.  By  substituting  Equations  E.ll  and  B.12 
into  Equations  B.1-B.9,  the  [C]  matrix  of  transversely 
isotropic  material  is  then  obtainel.  Equations  B.ll  and 
B.12  show  that  there  are  only  five  independent  material 
constants,  namely,  Ej^,  E,p,  and 

B.2  Formulations  of  Finite  Element  Method 

A  three-dimensional  elastic  body  bounded  by  surface  S 
is  loaded  by  concentrated  forces  p  on  nodes  and  surface 
traction  t  on  surface  SI.  Assume  that  the  body  force  is 
negligible  and  that  the  whole  system  is  adiabatic  and 
conservative.  According  to  the  principle  of  stationary 
potential  energy  (58],  the  equilibrium  is  ensured  if  the 
total  potentiaJ  energy  is  stationary  for  variation  of 
admissible  displacement,  l.e. 

6  (U  -  W„)  =  0  B.13 

where  U  is  the  strain  energy,  and  w^  is  the  work  done  by 
external  loads.  In  fact,  in  the  elastic  case,  potential 
energy  is  not  only  stationary  but  is  a  minimum  (59]. 

If  there  is  not  any  initial  stress  or  strain,  the 
strain  energy  U  and  the  work  done  by  the  external  load  are 


V 


U 


[  {c}  (o}  dV 


B.14 
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and 


W 


e 


{u}  dA 


where 


{e}  ~  [e^  -2  ^yz  ’*^xz 


A 

{o}  -  [Ojj  Oy  O2  ty2  1^2  "^xy^ 


{P}  -  [P^i  ?yi  P2i] 


{Ui}^  “  Uyj^  '^Zi^ 


T 

{u}  =  [Ujj  Uy  U^l 


T 

{t}  -  [t^  ty  t^] 


B.15 


B.16 


B.17 


B.18 


B.19 


B.20 


B.21 


in  which  is  the  nodal  displacement  in  x  direction  at 
the  ith  node,  and  is  the  coi.centrated  force  at  the  ith 
node  in  x  direction,  m  is  the  total  number  of  nodes  where 
concentrated  load  is  applied. 
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By  introducing  stiffness  matrix  [C]  and  the  shape 
function  matrix  [N],  i.e. 

{o}  =  IC]  {e}  B.22 

{u}  =  [N]  {q}  B.23 

strain  {e}  can  be  expressed  as  function  of  nodal 

displacement  {q}  by  using  strain-displacement  relatior: 

{c}  =  [L]  {q}  B.24 

Therefore,  Equations  B.14  and  B.15  can  be  rewritten  as 


1  T  T 

U  =  -  {ql  [Ll  [C]  (L]  (q}  dV 

V 


B.25 


m  T  r  T 

We  =  I  {P}  [NI  {q}  +  /  {t}  [N]  {q}  dA 

i-j  1  i  / 

®1 


B.26 


By  using  Equations  B.25  and  B.26,  a  system  of  3xn  equa¬ 
tions  is  obtained  from  Equation  B.13,  if  the  total  number 
of  nodes  is  n. 


[kl  {q}  =  {f} 


B.27 


155 


where 


T 

Ik],  ,  =  IL]  (Cl  (L)  dV 

3nx3n 


B.28 


^yl  ^zl  *  *  *  ^yn 

B.29 


{fl  =  I 

3nxl  i=j 


[Nl  + 


T 

{t}  [N]  dA 


B.30 


In  this  study,  the  interpolation  functions  of  linear 
and  parabolic  isoparametric  elements  (eight  nodes  and 
twenty  nodes,  respectively)  (55]  are  choccu  for  shape 
functions.  Each  element  has  the  same  thickness  as  the  ply 
thickness  of  fiber  composite.  Once  the  nodal  displf cement 
{q}  is  determined  by  solving  the  system  equations 
(Equation  B.27),  the  displacement  and  the  strain  fields 
can  be  obtained  by  Equations  B,23  and  B.24,  respectively. 
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